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Abstract 

A series of exact BPS solutions are found for single and double domain 
walls in M = 2 supersymmetric (SUSY) QED for finite gauge coupling 
constants. Vector fields are found to be massive, although it is localized 
on the wall. Massless modes can be assembled into a chiral scalar mul- 
tiplet of the preserved Af = 1 SUSY, after an appropriate gauge choice. 
The low-energy effective Lagrangian for the massless fields is obtained 
for the finite gauge coupling. The inter-wall force is found to be much 
stronger than the known infinite coupling case. The previously proposed 
expansion in inverse powers of the gauge coupling has pathological oscil- 
lations, and does not converge to the correct finite coupling result. 
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1 Introduction 



In recent years, much efforts have been devoted to the brane world scenario [l]-[3], where our 
world is supposed to be realized on a subspace such as a wall embedded in higher dimensional 
spacetime. On the other hand, supersymmetry (SUSY) has been one of the most fruitful ideas 
to build realistic unified models beyond the standard model [4]. SUSY helps to obtain domain 
walls and other stable solutions by requiring a part of SUSY to be preserved by the field config- 
uration. Then the configuration has minimum energy with the given boundary condition and is 
automatically a solution of field equations. Wall configurations can preserve half of SUSY charges 
[5]— [12] , which are called \ BPS states [13]. To obtain a wall with the four- dimensional world 
volume, we need to consider a fundamental theory in five or more spacetime dimensions. If such a 
theory are supersymmetric, it must have at least eight SUSY charges [14] . Simplest theories with 
eight SUSY consist of hypermultiplets. The nontrivial interactions of hypermultiplets require 
either nonlinearity of kinetic term (nonlinear sigma model) or gauge interactions [15]. The target 
spaces of such nonlinear sigma models are hyper-Kahler (HK) manifolds [16, 17] and can most 
conveniently be obtained by using vector multiplets as Lagrange multiplier fields [18, 19]. By 
introducing the Fayet-Iliopoulos (FI) term for the U(l) Lagrange multiplier vector multiplet and 
mass terms for hypermultiplets, one can obtain nontrivial potential terms with N discrete SUSY 
vacua for N hypermultiplets. Domain wall solutions have been obtained in these massive nonlin- 
ear sigma models with the target space T*CP N ~ l [20]-[26] and their generalizations [27]. One of 
the most convenient methods to obtain these nonlinear sigma models is to take a limit of infinite 
gauge coupling from the usual minimal gauge interactions of hypermultiplets, namely massive 
J\f = 2 SUSY QED. However, it has been known that the massive M = 2 SUSY QED itself has 
discrete SUSY vacua and may allow domain walls even with finite gauge coupling [28]-[32]. 

The massive nonlinear sigma models with multi-flavor have been found to give BPS multiple 
wall solutions containing a number of parameters, which represent the relative positions and 
phases of these multiple walls and are called moduli [28]. When these moduli parameters are 
promoted to functions of the world volume coordinates of the walls, one can obtain the low- 
energy effective Lagrangian of moduli fields, which are massless [33, 34, 22]. It has been noted 
that gauge field should acquire nontrivial configuration when these moduli depend on the world 
volume coordinates [30]. It has also been proposed that an expansion in inverse powers of gauge 
coupling may be useful to obtain the moduli dynamics at finite gauge coupling [28]. 

To fulfill the goals of the brane-world scenario, we need to localize particles of the standard 
model to "branes" , such as domain walls. It has been a notoriously difficult problem to obtain 
massless localized gauge field on a wall [2, 35, 36]. The general arguments [2, 35] and an explicit 
model [36] suggest that gauge fields may be localized on a wall, but will not be massless in the 
massive M = 2 SUSY QED with finite gauge coupling. Nevertheless, it is still interesting to study 
explicitly the role and the properties of dynamical gauge fields in the interacting system of the 
vector multiplet and hypermultiplets. In this respect, it should be useful to obtain exact solutions 
with finite gauge couplings [32] 1 . In the four-dimensional massive Af = 2 SUSY QED which 
can be obtained by a dimensional reduction of our model, it has been argued that a massless 
gauge field is obtained by dualizing the massless Nambu-Goldstone scalar [30] 2 . While this is an 

1 In Ref.[32], an exact solution of domain wall junction has been obtained for the massive Af = 2 SUSY QED 
with the finite gauge coupling for the first time. Therefore it automatically implies an exact solution of domain 
wall for finite gauge coupling. 

2 They are interested in an TV = 1 gauge theory obtained from J\f = 2 SU (2) gauge theory deformed by an 
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intriguing result, the gauge field as a dual of a compact scalar is an intrinsically three-dimensional 
peculiarity which is not valid to the more realistic situation such as our five- dimensional theory. 
In another paper, we show that massless localized gauge fields can be obtained by coupling our 
model to a tensor multiplet [37] . 

The purpose of our paper is to work out exact solutions of BPS single and multiple walls 
with finite gauge coupling in massive Af = 2 SUSY QED in five dimensions and to study the 
fluctuations on the domain wall background. We obtain a series of exact solutions for single as 
well as multiple walls provided gauge coupling has particular values relative to the ratio of the 
mass parameters and FI parameters. We also show that the gauge multiplet is massive for any 
finite gauge coupling, irrespective of the possible existence of exact solutions. This is consistent 
with the previous general arguments [2, 35] and explicit examples[36]. In our model, the massless 
Nambu-Goldstone scalar still exist, but are apparently unrelated to the dual gauge field. Instead, 
we obtain the Nambu-Goldstone scalars and the Nambu-Goldstone fermions which form a chiral 
scalar multiplet of the remaining Af = 1 SUSY, by choosing an appropriate gauge. Moreover 
we have only massive spectra for the vector field. These results are in contrast to the previous 
results of the massive Af = 2 SUSY U(l) gauge theory in four dimensions, where the Nambu- 
Goldstone scalar is identified as the dual of the three-dimensional gauge field [30] . As a result of 
exact solutions, we also obtain the low-energy effective Lagrangian for massless moduli fields for 
finite gauge couplings, by promoting the moduli parameters into four-dimensional massless fields 
[33, 34, 22]. In the process of carrying out the program, we observe that the gauge field should 
acquire nontrivial configuration with nonvanishing field strength around the wall, confirming the 
observation in Ref.[30]. Moreover, we find that this nontrivial gauge field configuration should 
be determined by means of the field equation of gauge field. We compute also some lower order 
terms of the previously proposed expansion in inverse powers of gauge coupling [28] , and compare 
the expansion with our exact result. We find that the expansion oscillates wildly and does not 
converge to the exact result in any smooth or uniform way. Therefore the expansion has a 
difficulty or at least subtlety in extracting physical quantities. 

In Sect. 2, we introduce our model of massive Af = 2 SUSY QED and obtain its SUSY vacua 
and BPS equations. In Sect. 3, we study the single wall in the case of two hypermultiplets 
(N = 2), and work out fluctuation spectra to show that there are massless Nambu-Goldstone 
modes, but that the vector field has only massive fluctuations. In Sect. 4, a series of exact 
solutions are constructed for N — 1 walls in the case of N hypermultiplets and obtain the low- 
energy effective Lagrangian for the massless fields. We also work out a few lower order terms of 
the expansion in inverse powers of gauge coupling and compare it to our exact solution. 

2 Massive N = 2 SUSY QED and BPS equations 

The simplest building block of Af = 2 SUSY theory in five dimensions is hypermultiplets which 
consist of SU(2) R doublets of complex scalar fields H tA , Dirac fields ip A and complex auxiliary 
field F t A , where i — 1,2 stands for SU(2) R doublet indices and A = 1, ■ ■ ■ , N stands for flavours. 
For simplicity, we assume that these N hypermultiplets have the same U(l) charge, say, unit 
charge. The U(l) vector multiplet consists of a gauge field Wm, a real scalar field S, SU(2) R 

adjoint scalar mass term. However, they observed the model reduces to our Af = 2 SUSY QED to the leading 
order of the mass deformation parameter. 
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doublet of gauginos A* and SU(2) R triplet of real auxiliary fields Y a , where M, N = 0, 1, • • • , 4 
denote space-time indices, and a = 1,2,3 denotes SU(2) R triplet index, respectively. In this 
work, we shall consider a model with the minimal kinetic term for hypermultiplets and vector 
multiplets. The Af = 2 SUSY allows only a few parameters in our model : the gauge coupling 
g, the mass of the A-th hypermultiplet m A) and the FI parameters ( a . The FI parameters are 
real and transforms as a triplet under SU(2)r. Then the bosonic part of our Lagrangian reads 
[28]-[32], [38] 

£boson = --\{F MN {W)) 2 + -^(<9 M £) 2 + (V M H)\ A (V M H iA ) - H\ A {Y> - m A fW A 
^9 ^9 

+ ^- 2 (Y a ) 2 - ( a Y a + Hj A (a a Y a Y^ A + F%F A , (2.1) 
^9 

where a sum over repeated indices is understood, F MN (W) = OmWn — B^Wm, covariant deriva- 
tive is defined as Dm = 9m + iWm, and our metric is t/mn = (+1, — 1, • • • , — 1)- 

To ensure the discreteness of SUSY vacua, we need to make all hypermultiplet masses non- 
degenerate. Without loss of generality, we assume the following ordering of the hypermultiplet 
mass parameters 

m A+1 < m A (2.2) 

for all A. When all hypermultiplet masses are nondegenerate, the symmetry of our model reduces 
to U(1) N . The diagonal U(l) is gauged by the vector multiplet Wm, and the remaining [/(l)^ -1 
is the global symmetry. 

The easiest way to find SUSY vacua is to explore the condition of vanishing vacuum energy. 
To facilitate the procedure, let us first write down equations of motion of auxiliary fields Y a and 

Ft 

Y a = g 2 [C-Hl A (a a y,W A ], (2.3) 
F A = 0. (2.4) 

After eliminating the auxiliary fields, we obtain the potential V pot 



£boson = -^(F MN (W)) 2 + — ^mX) 2 + iV M H)\ A {V M H lA ) - U pot , (2.5) 



V pot = 



_ 9 
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C-H\ A {a a y 3 W A +H} A (E-m A ) 2 H iA (2.6) 
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The vanishing vacuum energy is achieved by requiring 

hU^YjH^ = C, (2.7) 

for each a = 1,2,3 (the flavor index B is summed) and 

(E - m A )H iA = 0, (2.8) 

for each flavor A (the flavor index A is not summed). These SUSY vacuum conditions guarantee 
the full preservation of SUSY and can also be derived by requiring the SUSY transformation of 
fermions to vanish as we see immediately. 
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By making an SU(2)r transformation, we can always bring the FI parameters to the third 
direction without loss of generality 

C = (0,0,C), C>0. (2.9) 

In this choice, we find N discrete SUSY vacua (A = 1, • • • , N) explicitly by solving Eqs.(2.8) and 
(2.7) as 

S = m A , \H 1A \ 2 = (, H 2A = 0, 

H 1B = 0, H 2B = 0, (B^A). (2.10) 



Since fermions are assumed to vanish in the wall configuration, we need to examine only 
SUSY transformations of fermions to find a configuration preserving a part of SUSY. Gaugino 
X 1 and hyperino i\) A transforms as 3 

5 £ X = (I 7 ^F MW (PU)+7 M 9Ms)e l + ? (rV a )>, (2.11) 

5^ A = -iV2~U M V M H lA + i(S - m A )H iA ] e tj e 3 + V2F A e\ (2.12) 

To obtain a wall solution, we assume the configuration to depend only on the coordinate of one 
extra dimension, which we denote as y = x A . We also assume the four-dimensional Lorentz 
invariance in the world volume coordinates x^ = (x°, • • • , x 3 ), which implies 

F MN (W) = 0. (2.13) 

Since we are interested in the | BPS configuration, we require the above SUSY transformations 
(2.12) to vanish for half of the Grassman parameters specified by 

p +e x = o, P_e 2 = 0, (2.14) 

where P± = (1 ± 7s)/2 are the chiral projection operators. Finally we need to eliminate the 
auxiliary fields Y a and F A by their algebraic equations of motion (2.4) to make the BPS condition 
as first order differential equations for physical fields. Thus we obtain the | BPS equations for 
the massive M = 2 SUSY QED as 

d y Z = g 2 (C - Hl A H 1A + H\ A H 2A ) , (2.15) 

2g 2 H\ A H 1A = 2g 2 H\ A H 2A = 0, (2.16) 
V y H iA = (m A -J:)H iA , i = 1,2, A = 1,---,N. (2.17) 

One can easily see that the translation invariant vacuum requires the vanishing of the left-hand 
side of Eqs.(2.15) and (2.17), which implies the same condition as the full preservation of SUSY 
in Eqs.(2.8) and (2.7). 

We are interested in the solution of these BPS equations which interpolate two different vacua 
in Eq.(2.10). Since the BPS equation is a first order differential equation, the boundary condition 
dictates that H 2A should vanish identically 

H 2A (y) = 0. (2.18) 



3 Our gamma matrices are 4x4 matrices and are defined as : {j M , 7^} = 2rj MN , ^ MN = i [y M , j N ] = j\- m j n ^ , 
7 5 = Z7°7 1 7 2 7 3 = — ij 4 . 
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The energy density of the BPS solution can be found by making a Bogomolny completion of 
the energy functional 



Tw = II dy iM dyj: " g2{c ~ \ HlA ^) 2+ \ v y HlA - ^ - m 1A \ 2 } 

+ [CS + (m A - Y)\H lM ~ 



y=+oo 
y=-oo 



(2.19) 



yielding the BPS wall tension for solutions interpolating the A-th vacuum and the _B-th vacuum 

T w = ((m A -m B ), (2.20) 

assuming A> B. The structure of the above BPS equations shows that only those hypermultiplet 
scalars H 1C ^ with A > C > B have nonvanishing values besides the vector multiplet scalar 
£ [22], [27]. Defining I = A — B, we call such a BPS configuration as I wall solutions, since it 
represents / separate walls at least when these walls are sufficiently far apart, as we shall see in 
the Sect. 4. 



3 Single wall BPS solutions 

3.1 Exact solutions for single wall 

Since single wall interpolates two adjacent SUSY vacua, we shall assume N = 2 (two hypermulti- 
plets), without loss of generality. For simplicity, we take mass parameters of the hypermultiplets 

as 

mi = —m>2 = m - (3.1) 

The boundary conditions for BPS equations are given by 

£(— oo) = —to, S(oo) = to, (3.2) 
if n (-°o)=0, if^oo) = v/C, 

H l2 (-oo) = y/C, H 12 (oo) = 0. (3.3) 



The set of the BPS equations (2.15)-(2.17) with the above boundary conditions (3.2)-(3.3) 
are known to be solved exactly for infinite gauge coupling [22], [27]. Recently another exact 
solution has been found for finite gauge coupling, provided the gauge coupling g satisfies the 
following relation with the ratio of the FI parameter ( and the mass difference squared (2to) 2 
[32] 

g 2 ( = 2m 2 . (3.4) 

By generalizing this exact solution for finite gauge coupling, we find that there are exact solutions 
for the following values of gauge coupling g 

it - 8m2 <\ 
9 C = -p-, (3-5) 
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with appropriate integers k. We shall denote the single wall solution for k as Sk{m). The integer k 
indicates the value of the gauge coupling g through the relation (3.5), and the mass parameter m 
indicates the tension T w of the single wall in unit of the FI parameter as T w = 2(m. The infinite 
gauge coupling corresponds to the case k = 0, and the previously obtained finite gauge coupling 
case to k = 2 [32]. We have explicitly obtained exact solutions in the case of k — 0, 2, 3, 4. 

In this section, we will describe the exact solution 5*2 (m) for k = 2 as the simplest example, 
which plays an important role in our paper on massless localized vector field [37] : 

= mtanhrm/, (3-6) 

777 

H ll {y) = — ^e +my sech my, (3.7) 
V2g 

IT) 

H 12 (y) = -— e - my sech my. (3.8) 
V2# 

Let us recall that we can introduce two arbitrary integration constants in these solutions. We can 
have y — > y — yo corresponding to the spontaneously broken translation invariance, and we can 
also have the multiplication by a phase e~ ia ( e ia ) for hypermultiplets if 11 (if 12 ) corresponding 
to the spontaneously broken global U(l) invariance. They constitute two moduli of the solution. 
The yo has a physical meaning of the position of the wall. The vector multiplet scalar E and 
the hypermultiplet scalars if 11 and if 12 of the solution 5*2 (m) are illustrated as a function of the 
coordinate y in the extra dimension in Fig.l. We will present other exact solutions of single wall 
after introducing a slightly different notation suitable also for multiple walls. 



S/m V2gH 1A /m 




a) b) 

Figure 1: The exact BPS wall solution S^ra) for finite gauge coupling g = \/2m/^/C, and the 
tension T w = 2m(. a) Scalar field E(?/) of vector multiplet in Eq.(3.6) divided by the mass 
parameter m as a function of the coordinate y times m. b) Hypermultiplet scalar field H u (y) 
in Eq.(3.7) as a function of the coordinate y times m (solid line), and H 12 (y) in Eq.(3.8) as a 
function of the coordinate y times m (dotted line). 



3.2 Fluctuation around the BPS single wall background 

Let us now turn our attention to spectra of fluctuations around the BPS single wall background. 
These fluctuation fields can be expanded into modes defined on the background. Among various 
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modes, we are particularly interested in two kinds of modes : massless modes (zero modes) 
and fluctuation of vector fields. The former modes play essential role in discussing dynamics 
on the background at low-energies. The latter is to examine the role played by the dynamical 
vector field at finite coupling, rather than the role as the Lagrange multiplier field in the case of 
corresponding nonlinear sigma models. 

Zero modes usually arise from several reasons. The first one is the Nambu-Goldstone modes 
as a result of the spontaneous breaking of continuous global symmetry. The second possibility is 
the result of remaining symmetry, such as the conserved SUSY [39] . The third possibility is the 
result of index theorems, such as fermion zero modes around instantons and other solitons [40]. 

To illustrate the method, we shall work on the fluctuations on the background of the simplest 
exact solution S^m) in Eqs.(3.6)-(3.8). In the case of a single wall, such as S^m), the only 
spontaneously broken global internal symmetry is the U(l) symmetry which rotates the phase 
of hypermultiplets oppositely 

H il _^ e ~i* H i^ H i2 _, e ia#*2 ( 3 9 ) 

There are two bosonic global symmetry which are spontaneously broken : translation and the 
global U(l) in Eq.(3.9). The mode function of the Nambu-Goldstone boson corresponding to the 
spontaneously broken translation is given by differentiating the background field configuration 
with respect to the coordinate in extra dimension y. We shall denote the corresponding four- 
dimensional effective field as Re0 o (^)- Similarly the Nambu-Goldstone boson associated with the 
global U(l) symmetry is given by acting an infinitesimal U(l) transformation. We denote the 
corresponding four-dimensional effective field as lmcf) (x). By writing out only zero modes, we 
thus obtain scalar fields of vector multiplet E(x, y) and that of hypermultiplets H 1A (x, y),A— 1,2 

Z(x,y) = (£(y)> -^(£(y))Re0o(z), (3.10) 
H u (x,y) = (H u (y)) - d y (H n (y))ReM^ - im(H 11 {y))lM {x) , (3.11) 
H 12 (x,y) = (H l2 (y)) - d y (H l2 (y))ReMx) + im{H 12 {y))l^ct> {x) , (3.12) 

where we multiplied the factor m in front of the U(l) Nambu-Goldstone mode wave function for 
later convenience, and the (• • •) denotes the background field configuration. These configurations 
together with vanishing gauge field Wm = are solutions of the BPS equation if Re^o and Im0 o 
are constants (moduli) and hence automatically satisfy the field equations. However, we are 
now interested in the case where these Re0 o and Im0 o are four-dimensional fields, which depend 
on the coordinate x M of the world volume. Then these field configurations are no longer BPS 
configurations, nor satisfy the field equations. Since four-dimensional effective fields are defined 
as solutions of the field equations, we should demand these field configurations to satisfy the 
linearized field equations. Since Wm must vanish for constant <po, it should be proportional to 
d^<f)o(x). Therefore the Lorentz invariance in the four- dimensional world volume dictates 

W 4 (x,y) = 0, (3.13) 

and 

W li (x,y)=w(y)d l Jm<l>o(x)- ( 3 - 14 ) 
with some function w(y) to the linearized order which we are interested in. Then the linearized 
field equations for the vector multiplet scalar E, hypermultiplet scalars H lA are satisfied by 

d^Mx) = 0, (3.15) 
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which just means that these Nambu-Goldstone fields are massless. 

The remaining linearized field equation for the gauge field Wm reads 

d*PWp - d 2 W, - d^drW) + d,d y W A = -2g 2 {(H 1A )d,lmh 1A + \(H 1A )\ 2 W,) , (3.16) 

dvVWt - d y {d v W v ) = 2g 2 (d y (H 1A )lmh 1A - (H 1A )d y lmh 1A - \(H 1A )\ 2 W 4 ) ■ (3.17) 

The solutions for Wm in Eqs.(3.13) and (3.14) are enough to satisfy the field equation for W4 in 
Eq.(3.17). The remaining field equation (3.16) for determines the y dependence w(y) of the 
gauge field W^x^y) in Eq.(3.14) as 

W^x,y) = &{y))d l M4> Q {x). (3.18) 

This solution shows that the corresponding field strength is nonvanishing around the wall 

F^(W)(x,y) = -d y (E(y))d l JmMx). (3.19) 

In a similar model in four- dimensions, it has been observed that the gauge field must be nontrivial 
and field strength must be nonzero around the wall if the moduli depends on world volume 
coordinates [30, 28]. Our result in five dimensions are in agreement with their observation. Our 
new point is perhaps that the necessary gauge field configuration can be determined explicitly 
by using the linearized field equations for gauge fields. 

SUSY is another global symmetry which is spontaneously broken. We can obtain the Nambu- 
Goldstone fermion corresponding to the broken SUSY by evaluating the SUSY transformations 
of fermions on the background. One set of S'6 r (2)-Majorana spinor A 1 , % — 1, 2 in five dimensions 
can be decomposed into two Majorana spinors A + , A_ in four dimensions as 

A + = P+A 1 + P_A 2 , A_ = -P_X 1 + P+A 2 . (3.20) 

The SUSY transformations in Eqs.(2.11) and (2.12) can be evaluated by using the BPS equations 
(2.15)-(2.17) to yield for the gauginos 

5 £ \+ = i [d y {X) +g 2 {(- \(H lA )\ 2 )] (P + e + - P_e + ) = 2id y {T) {P + e + - P_e + ) 

5 £ \- = i [d y (T)-g 2 (C- \(H 1A )\ 2 )] (P_ £ _ +P +£ _) =0. (3.21) 

Since e_ is preserved, it drops out from these transformations. If we evaluate explicitly for our 
exact solution S 2 (m) in Eqs.(3.6)-(3.8), we obtain 



= —2 — ( P + £ + ~ P - e +) > = °- (3-22) 



2im 2 
cosh 2 my 

Similarly the hyperino component of the Nambu-Goldstone fermion wave function is given by 

5 £ ^ A = V2 [d y (H 1A ) + ((E) - m A ) (H 1A )} P + e_ + V2 [-d y (H 1A ) + ((E) - m A ) (H 1A )} P_e + 
= -2V2d y (H 1A }P_E + . (3.23) 

For the exact solution S 2 (m) in Eqs.(3.6)-(3.8), it becomes explicitly as 



The above transformation laws (3.21) and (3.23) show that — 4e + is proportional to the zero 
momentum component of the Nambu-Goldstone fermion. Therefore we can define the Nambu- 
Goldstone fermion field Xoi x ) °f the four-dimensional effective low-energy theory on the world 
volume of the wall as 

K(x,y) = - l -d y (X(y)) l5Xo (x), \-(x,y) = 0, (3.25) 

^ A (x,y) = -±=d y (H 1A (y))P_ Xo (x), (3.26) 

suppressing to write massive modes. We see explicitly that the Nambu-Goldstone fermion appears 
only with the left-handed chirality in the hypermultiplets as is usually dictated by index theorems 
for fermions localized on domain walls [40]. We can also verify that the Nambu-Goldstone fermion 
satisfies the linearized equations of motion with a vanishing mass eigenvalue, by using the BPS 
equations (2.15)-(2.17). 

Now let us consider the requirement of the preserved symmetry, especially SUSY. The SUSY 
transformation property under the preserved SUSY e_ specified by Eq.(2.14) is given by 

Se.W^ = -ie_ 7 „A_, (3.27) 

5 £ _W A = e_A+, (3.28) 

5 £ _S = -i£_ 75 A+, (3.29) 

5 £ _H 1A = -V2e-P-ip A , (3.30) 

5 £ _H 2A = -V2e_P + *lj A . (3.31) 

This transformation property does not fit well with the above massless particles of Nambu- 
Goldstone modes in Eqs(3.10)-(3.14), (3.25), and (3.26), especially for the vector multiplet. 
However, we should remember that we have a freedom to make gauge transformations to make 
the above massless particles in conformity with the SUSY transformation properties. Let us 
perform the following gauge transformation, which is proportional to the fluctuation field Im0 o 

W 4 — > W 4 - d„ «£)Im0o) , 
W^^^-^((£)Im0 o ), 
H 1A — (H 1A ) + i ((£)Im0o) (H 1A ). (3.32) 

Then the hypermultiplet scalars become 

H u (x, y) = (H u (y)) - d y (H l \y))ReUx) -i(m- (E(y))) (H 11 (y))lm0 o (x) 

= (H u (y)) - d y (H u (y)) (Re0 o (x) + Om^x)) , (3.33) 

H 12 (x, y) = (H 12 (y)) - d y (H 12 (y))ReMx) -i(-m- (E(y))) (H 12 (j))H(i) 

= (H 12 (y)) - d y (H 12 (y)) (Re0 o (x) + Tm0 o (^)) ■ (3.34) 

Therefore we find that the two Nambu-Goldstone bosons corresponding to translation and U(l) 
global rotation forms a complex scalar <po(x) of a chiral scalar multiplet of the preserved four 
SUSY e_. At the same time, the four- dimensional gauge fields vanishes, and vector multiplet 
scalar S and the fourth (extra dimension) component of vector field W 4 fit into the complex 
scalar <po(x) 

X(x,y)+iW 4 (x,y) = (S(y))-^(S(y))(Re0 o (x)+2lm0 o (x)), (3.35) 
W li (x,y) = 0. (3.36) 
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These mode expansions together with the fermionic ones in Eqs.(3.25) and (3.26) are now consis- 
tent with the preserved SUSY which defines the four-dimensional M = 1 SUSY transformations 
for the massless fields in the effective field theory 

6 e _Mx)=e-P-Xo(x). (3.37) 

Thus we find that all the Nambu-Goldstone bosons and fermions together form a chiral scalar 
multiplet (0o (#)) Xo( x )) of the preserved SUSY. Moreover we do not find any index theorem to 
force us additional zero modes. Therefore we do not expect any more zero modes. 

Let us now turn to the issue of mass spectra, especially for vector fields. The linearized 
equation of motion for vector field and PU 4 are given in Eqs.(3.16) and (3.17). By exploiting 
the U(l) gauge invariance, we can always choose a gauge where W±{x,y) = identically. We can 
also decompose vector field into the transverse part and longitudinal part w 

W li (x,y)=W li (x,y) + d li w(x,y)> &% = Q. (3.38) 

The transverse part of the vector field 4 is decoupled from the rest of the fluctuations because 
of the four-dimensional Lorentz invarince and gives the following mode decomposison with the 
mode function a n (y) for the mass eigenvalue m n of the four- dimensional effective field wj?\x) 

W fM (x,y) = J2^n(y)W^(x), (3.39) 

n 

-d 2 a n (y) + V(y)a n (y) = m 2 n a n (y). (3.40) 

where the potential V(y) is defined by 

V(y) = 2g 2 \{H 1A )\ 2 : (3.41) 

and is illustrated in Fig. 2 for the background ^(m) in Eqs.(3.6)- (3.8), and £3(771) and £4(777), 
whose exact solutions are given explicitly in Sect. 4.1. These potentials have a dip near the 



wall and approaches 2g 2 ( at y — > ±oo. Therefore all the vector fluctuations become infinitely 
heavy as g — > oo. Potentials are always positive definite, but have attractive forces due to a dip 
around the wall for any finite gauge coupling. It is likely that there may be discrete low mass 
states before reaching a continuum starting from m 2 = 2g 2 (. The potential becomes infinite for 
the infinite gauge coupling. It is clear that the eigenvalue can never vanish : m n > 0. It is also 
easy to see that m n are of the order of m for, say, £2(777). This result is consistent with the 
general argument that the gauge field will obtain a mass of the order of the inverse width of the 
wall. Although the low-lying vector modes have presumably discrete spectra and are localized 
around the wall, all these modes are massive. To obtain a massless vector field, we can add 
tensor multiplet [37]. 

To find out the remaining field equations, we define fluctuations of field as E = (E) + s, 
H A% = {H Al ) + h Al . The linearized field equations are given by 

Q = d M d M s + 2g 2 {\H M \ 2 )s + 2fl 2 /4((£ - m A )H M ) + 2< 7 2 <if* A (£ - m A ))h A \ (3.42) 

4 This component corresponds to the polazation states proportional to the polazation vector e^,i — 1,2,3 
defined by e^p M = in momentum space, representing a genuine four-dimensional vector field. 
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V(x) 




Figure 2: The potential V(x) of Eq.(3.41) for ^(m) (solid line), Ss(m) (dotted line), and £4(771) 
(dashed line) in the eigenvalue equation for the vector fluctuation (3.40) as a function of extra 
dimension coordinate y. The mass parameter is taken to be m = 1. 

= d M d M h Ai + i(d M W M )(H Ai ) + 2iW M d M (H M ) - g 2 (( - (\H Bj \ 2 ))h Ai 

+ g 2 (H M ) ((rf Bj )h B i + (H* B )h* Bj ) + 2«E> - m A )(H Ai )s + ((E) - m A fh Ai . (3.43) 

These linearized field equations can be decomposed into three sets of coupled equations. The 
first set consists of the fluctuations of vector multiplet scalar E, and those of the real part 
of the hypermultiplet scalars ReH 1A , A = 1,2. The second set consists of imaginary part 
of the hypermultiplet scalars ImH 1A , A = 1, 2 and the longitudinal part of the vector field 
w in Eq.(3.38). The third set consists of the lower components of the hypermultiplet scalars 
H 2A , A = 1,2. Among various massive modes, we find the following tower of massive modes as 
a solution of the first set of coupled equations 

s = J2( d y a n(y))U^), h Al = - J2( HA1 (y)) a n(y)Mx), (3.44) 

n n 

with h = Wm = 0. This effective field <p n (x) has the identical mass squared vr^ n as the effective 
field wjr\x) of the transverse vector in Eq.(3.40). Moreover the mode function is also given by 
the same wave function a n (y) as the transverse vector. Therefore it is likely that these effective 
fields are related by the preserved symmetry, in particular the Af = 1 SUSY. In fact, a massive 
vector multiplet of Af = 1 SUSY should contain a scalar particle. 

4 BPS solution for Multiple walls 

4.1 Exact solutions with finite gauge couplings 

It is convenient to introduce a complex function tp(y) to solve the BPS equation for hypermulti- 
plets (2.17) [28] 

N-2 

H 1A {y) = VC™p(-^(y)+rn A (y-y ) + Y,®A a r a ), (4.1) 

a=l 

where the complex parameters r a ,a = 1, ■ ■ ■ , N— 2 are collective coordinates arising as integration 
constants. Since two complex parameters among the integration constants can be absorbed by 
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a shift of tp and y , we can choose the N x (N — 2) fixed real matrix a A a to be of rank N — 2. 
The BPS equation for hypermultiplets is equivalent to the following equation for ip(y) 

d y ip = i: + iW A . (4.2) 

Because of the vanishing field strength (2.13), the vector field W4 = Imdyip is a pure gauge. 
However, when we consider dynamics of domain walls, this term will play an important role. 
We have the following moduli parameters in the solution (4.2), 

y = Y + i9 , (4.3) 
r a = R a + i6 a , (4.4) 

with one and N — 2 complex dimensions, respectively. The Y and R a are related to the center 
of mass and relative positions of the a-th domain wall, respectively, and 9 and 9 a to the overall 
phase and relative phases of the a-th wall, respectively. Using the variable ip, the BPS equation 
for vector multiplet (2.15) becomes 

N N-2 

—d 2 y Re(ij) = 1 - ^exp(-2Re(^) + 2m A (y - Y ) + 2 ^ a A a R^) 

^9 a=1 a=l 

= 1 - exp (-2Re(V>) + 2W) . (4.5) 
Here, the explicit dependence on y in Eq.(4.5) can be assembled by defining a function W(y) as 

N N-2 

W = log ex P ( 2m A(y -Y )+2Y, u A a R a ) ■ (4.6) 

A=l a=l 

To ensure that d y T, = at y — > ±00, we should impose boundary conditions for the above 
equation as 

N-2 

Re(V>) — > m 1 (y - Y ) + ^ afRa, y -> +00, (4.7) 

a=l 

N-2 

Re(^) — ► mjv(y - y ) + Yl a N aR a, V -> -00. (4.8) 

a=l 

To define an appropriate variable for the position of each wall, let us consider a configuration 
where only two adjacent exponential terms in the sum of Eq.(4.5) are large and the others are 
negligible, then the function W has a profile for a single wall. Therefore it is natural to define 
the position y A of the A-th wall as 

(N-2 \ / N-2 

2m A (y A - Y ) + 2 Y ®A a R a J = exp I 2m A+l (y A - Y ) + 2 a^i"^ 
a=l / \ a=l 

- ^ = y - E - l2K+ia ~^ a)jRa . (4.9) 

m A+1 - m A 

This moduli parametrization y A has an intuitive meaning of the position of the A-th wall and 
yA — Va+i corresponds to a distance between two walls, at least when the distance is large. 
Relations between y A and the relative positions R a defined in Eq.(4.4) are obtained by 

N-2 A-l 

- Y,( a A a - a ± a )R a = (rn B+1 - m B ) (y B -Y ), A = 2, • • • , N. (4.10) 

o=l B=l 
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By choosing A = N, we find 

1 



m,N — mi 



which becomes the center of mass coordinate when a a N 



N-l N-2 

S ^{m B+1 -m B )y B + s ^(a a N -a a l )R a , (4.11) 

B=l a=l 

a\. Because of the translational 
invariance, Re(^) should be a function of real variables y — Y and R a . From now on, we will 
take Yq — unless otherwise stated. 



Now we will present a series of exact solutions for single and double walls for finite values 
of gauge coupling. To compare single and multi-wall solutions, let us assign the following mass 
parameters for hypermultiplets for double walls (N — 3) : 

m A = {m,0, -m). (4.12) 

This convention is intended to make the total energy density of the double wall to be identical 
to the single wall, so that the double wall situation can be most naturally compared to the single 
wall situation with the same energy density (tension) in the coincident limit of two walls, since 



the mass parameters of the single wall is assigned to be m J 



to, 



-m) in Eq.(3.1), and the total 



energy density is just given by the difference of the two extreme masses as given in Eq.(2.20). 
For the moduli parameter, we use y\ — y<i = R 1 = R and choose the a A x as 

a A x = (0, to/2, 0), 2a A a R a = 2a A 1 Ri = (0, mR, 0), (4.13) 

since the rank of the matrix a A a is iV — 2. This relative distance appears only in multiple wall, 
but not in the single wall. 



The function W = Wangle for the single wall case becomes 

^ single ^ilog(e 2 ^ + e- 2 ^). 



Similarly the function for the double wall case W = WWbie is found to be 

double = \ log(e 2 ^ + e" 2 "^ + e mR ). 

The solvable cases of finite gauge coupling are found to be 

8to 2 



k 2 



k = 0,2,3,4, 



(4.14) 



(4.15) 



(4.16) 



where the mass parameter to is defined for single wall case in Eq.(3.1), and double wall case 
in Eq.(4.12). We will denote the single wall solution as Sk{m) and the double wall solution as 
Dk(m) with the coupling defined by (4.16) with k and the mass parameter to. 

Let us list the exact solutions that we are able to obtain 

Jlog(e 2m?/ + e- 2 "^), 



S (m) 


: Re^ 


S 2 (m) 


: Re^ 


S 3 (m) 


: Re^ 


S±{m) 


: Re^ 


D {m) 


: Re^ 


D 4 (m) 


: Re^ 



lo g (e mv + e- mv ), 
-log(ef m2/ + e~^ my ) 



\og(e my + e~ my 
1 
" 2 

-. log(e my + e~ my + \/6 + e 



^-\og(e 2my + e~ 2my + e mR ), 



mR 



)■ 



(4.17) 
(4.18) 
(4.19) 
(4.20) 
(4.21) 
(4.22) 
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It is interesting to observe that the double wall solution D (m) in Eq.(4.21) in the case of infinite 



mR 1 Q 




8 
6 




4 




(4.26), 




^/(4.24) 



-5 -2.5 2.5 5 7.5 10 



mR 



Figure 3: Relation between R defined in Eq.(4.4) and R. The solid line represents the relation de- 
fined in Eq.(4.26) for the finite coupling solution D±{m) in Eq.(4.22). The dashed line represents 
the relation defined in Eq.(4.24) for the infinite coupling solution Do(m) in Eq.(4.21). 



coupling can be rewritten as a superposition of two single walls placed apart by a distance R 

Re^(y) = - \og(e m(y -^ + e"" 1 ^ )) + - log(e m(j/+ f > + e - m{y+ ^). (4.23) 

This parameter R can be regarded as another choice of a moduli parameter and is related to R 
defined in Eq.(4.4) 



e mR + e - m R =e mR ^ m R = \og(- ±11 ( 4 .24) 



which is illustrated in Fig.3. Unfortunately the new choice of the moduli parameter R becomes 
pure imaginary and loses an intuitive meaning as the distance between the two walls when 
R < log2/m. This situation has some similarity to a moduli parameter for a model of BPS 
double wall in an M = 1 SUSY Wess-Zumino model in four- dimensions [41]. We can also rewrite 
the double wall solution D 4(171) in Eq.(4.22) in the case of the finite coupling as a superposition 
of two single walls placed apart by a distance R 

ReiP(y) = log(ef iy -^ + e^fo"! )) + log(e^ y+ ^ + e^^). (4.25) 
The new moduli parameter R is related to R in this case as 



e — +e _ — = V6 + e mR mR = 2 log^- — J. (4.26) 

In this case of the finite coupling D 4 (m) in Eq.(4.22), the parameter R takes only positive real 
values, whereas R takes positive as well as negative values as illustrated in Fig.3. In both cases of 
infinite and finite coupling, both R and R have an intuitive meaning of relative distance between 
the wall, as long as the distance is large : R — > R, for R — > oo. Therefore R in this case gives an 
intuitively nicer parametrization for the relative distance between the two walls. 
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Figure 4: Comparison of the vector multiplet scalar E as a function of my for exact solutions 
with various gauge couplings, a) Single wall solutions So (m) (dotted line), S2 (m) (solid line), 
S 3 (m) (short dashed line), and S 4(171) (dashed line), b) Double wall solutions D (m) (dotted 
line), and D 4(171) (solid line). 

It is interesting to examine how these double wall solutions behave, in the limit of coincident 
walls and in the limit of asymptotically far apart walls. By smoothly changing the moduli 
parameter R, we find the following limiting behaviors which can be expressed symbolically as 

S (m) R ^ D (m) So (|) S (|) , (4.27) 

S4(m) R ^° D 4 (m) R ^S 2 (^)®S 2 (^). (4.28) 

These limiting behaviors apply not only for the ip(y), but also for all the physical quantities. 
Note that the sum of the numbers k is preserved, corresponding to the fact that the total tension 
is a conserved topological invariant for a given boundary condition. We can find E, H 1A and the 
potential V corresponding to these exact single and double wall solutions (4.17)-(4.22) by using 
(4.1), (4.2) and (3.41). We illustrate and compare various single and double wall solutions for 
scalar E(y) of vector multiplet in Fig. 4, for hypermultiplet scalars H 1A (y) for A = 1,2 (single 
wall) and A = 1, 2, 3 (double wall) in Fig. 5 and for potential V(y) in Fig. 6. 




-6 -4 -2 2 4 6 -10 -5 5 10 



a) b) 
Figure 5: Comparison of the hypermultiplet scalars H lA as a function of my for exact solutions 
with various gauge couplings, a) Single wall solutions So(m) (dotted line), S 2 (m) (solid line), 
5*3 (m) (short dashed line), and S 4 (m) (dashed line), b) Double wall solutions D (m) (dotted 
line), and D 4 (m) (solid line). 
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a) b) 
Figure 6: Comparison of potential V = normalized by the asymptotic value (y — > ±00) as a 
function of my for exact solutions with various gauge couplings, a) Single wall solutions So(m) 
(dotted line), S^m) (solid line), Ss(m) (short dashed line), and S 4,(70) (dashed line), b) Double 
wall solutions D (m) (dotted line), and D±(m) (solid line). 



4.2 Moduli Dynamics of Two Domain Walls 

The low-energy effective Lagrangian for massless field has been studied systematically, especially 
in the case of infinite gauge coupling [28]. It has been proposed to use an expansion in inverse 
powers of gauge coupling, since the infinite gauge coupling case was the only explicitly known 
solution at that time. We would like to derive the low-energy effective Lagrangian for massless 
field in the case of finite gauge coupling, using our exact solutions. We will compare it to the 
previous result of infinite gauge coupling, and show an instability of an expansion in inverse 
powers of gauge coupling. 

Since the overall position and phase become just free massless fields [28], we shall concentrate 
on the relative distance R a and relative phase 9 a of multiple walls defined in Eq.(4.4). These 
moduli parameters represent positions in y and in internal space of U(l) of multiple walls. 
We have obtained mode functions of massless Nambu-Goldstone fields by differentiating the 
background in terms of these moduli parameters in Sect. 3. 2. Here we are interested in not only 
the mass spectrum but also the entire effective Lagrangian at low energies. To find out the 
low-energy effective Lagrangian, a more systematic method by Manton [33] is useful. In this 
method, we need to promote the parameters of the solution (moduli) to be functions of the 
world volume coordinates namely four-dimensional fields. Assuming the variation in the 
world volume coordinates x^ to be weak, we can obtain all the nonlinear dynamics containing 
smallest number of derivatives (two derivatives). This well-defined procedure allows us to obtain 
the nonlinear kinetic term of the massless fields. Since we are considering the massless Nambu- 
Goldstone fields, there should be no additional potential terms without derivatives. We can keep 
only t = x° dependence of these fields, since the Lorentz invariance in the world volume allows to 
recover the entire dependence on x^ 1 [28] . Effective Lagrangian C e s of moduli fields is found to be 
given by integrating over the Lagrangian of the fundamental theory evaluated by the background 
field with the t dependent moduli r a (t) = R a (t) + i9 a (t) [28] 

poo n 

£ cff = / dyJ2 l# 1A | 2 (^ + m A y - a a A r a )(m A yl - a A f\). (4.29) 

J -°° A=l 

Let us take two wall case (N = 3) for concreteness. We discard terms for the overall position 
and phase, and omit the subscript a, since there is only one set of relative moduli R, 9 now. We 
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will be interested in low-energy effective Lagrangian for the relative distance R and the relative 
phase 9 of two walls. To obtain the dependence on the relative phase 9, we need to find the 
imaginary part of i/j(y) besides the real part that we have determined. The BPS equation (4.2) 
shows that the imaginary part of ip(y) is given by integrating the extra dimension component 
of the gauge field : d y \m.il)(y) = W 4 (y). As shown in the case of the single wall in Sect. 3. 2, we 
have, near the walls, a nontrivial field strength F m (W) ^ proportional to the derivative of the 
moduli fields 9. Then the extra dimension component of the vector potential is determined by 
the field equation for the vector potential [28] . The most convenient gauge to make the remaining 
SUSY manifest is W M = and W 4 ^ 0. 

We find from the equation of motion that the imaginary part (combined with the real part) 
satisfies 

g = VC ^_,^±l!)) e — (4 .30) 

Because of the boundary condition ip — at y — ±oo, the solutions of the real and imaginary 
parts are proportional to R and 9 with the same coefficient which is a function of R only 

implying W 4 = d y -^Keip{y,R)9. The moduli field r = R + i9 is a complex scalar of a chiral 
scalar field, taking values [28] R e R and 9 G [0,47r/ra). We see that the M = 1 SUSY in 
four dimensions is now manifest. Moreover, the Kahler metric JC rr * of the low-energy effective 
Lagrangian C e s is given in terms of F(R) which is a function of R only 

C cS = (i? 2 + fl 2 )/C rr *(r,r*), (4.32) 
]C„,(r,r*) = -m(F{mR) = -m(F ( m, 1 —^— J . (4.33) 



4 * v ' 4 * V 2 

The function F(R) is most conveniently evaluated by taking 9(t) =0 and by writing hypermul- 
tiplets H 1A in terms of the function tp(y) in Eq.(4.1) and using the choice a a in Eq.(4.13) 

-m(R 2 F(mR)= / dy ^ \H 1A \ 2 (-Reipa A R + a 2 A R 2 ) 

A=l 

1 r\ r°° 

= t-mCR 2 ^ / dye~ 2 ^ +mR . (4.34) 
4 dR J.^ 

The Kahler metric /C rr * gives equations of motion for the relative distance R and for the 
relative phase 9 

R = - l -m{\ogF)\R 2 - 9 2 ), (4.35) 

6 = -m(\ogF)'R9, (4.36) 

respectively, where ' means d/dz, z = mR. These equations of motion are real and imaginary 
parts of a single equation due to the complex structure of the J\f = 1 SUSY : 

f = -\m{\og F)'r 2 , r = R + i9 (4.37) 
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We observe that the above equations of motion admit that 9 = is always a consistent solution. 
Therefore we can discuss the relative motion of the double wall with fixed relative phase consis- 
tently. On the other hand, in order to have a fixed R, the above equations of motion is consistent 
only if 9 is also constant. Therefore we shall consider the relative motion assuming a constant 
relative phase 9 in the following. 

Inserting our exact wall solution D^{m) with the finite gauge coupling in Eq.(4.22) into 
Eq.(4.34), we obtain the corresponding Kahler metric F(z = mR) for the finite gauge coupling 



dt- 

oo \ 6 



e* + e - * + 



P + e~* + V6 + e z ) 3 
e z - 4 + 



log 



y/6 + & + y/2T 



(4.38) 



(2 + e*)H v / ( 2 + e2 )( 6 + e2 ) 

For the D (m) solution at the infinite coupling, the Kahler metric F has been obtained as [28] 



-10 




Figure 7: Comparison of Kahler metrics F(z) for the infinite gauge coupling (dashed line) and 
F{z)\ for the finite gauge coupling (solid line), as functions of z — mR. 



~2 



f 

J — c 



dt- 



e* + e~* 



e* + e - * + e z ) 2 
4 



^2z 



e^ + 



ye 



2z 



log 



e 2 + y/e 



2z 



(4.39) 



The metric is real and positive for the entire values of — oo < z < oo, in spite of the apparent 
singularity at e 2z = 4 [28]. These Kahler metrics are illustrated and compared in Fig. 7. We see 
that there is a significant difference between them quantitatively, although general features are 
similar. 

The physical significance of the quantitative difference can perhaps most easily appreciated 
by comparing the strengths of forces acting between the walls in the case of the finite coupling to 
that of the infinite coupling. In Fig.8, we illustrate and compare the coefficients F'/F of force in 
Eqs.(4.35)-(4.37) for the infinite and finite coupling cases. We observe that the strength of the 
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force for the finite coupling case is much larger than the infinite coupling case for large relative 
distance R. 



F(z)'/F(z) 

if 




Figure 8: Comparison of the ratio F' Q / F (dashed line) and F'JF± (solid line). 



To understand the quantitative difference of the strength of the forces at large values of R, 
let us evaluate the Kahler metric F at asymptotic values of R. In the region of R — > — oo, both 
target space metrics F(z) are flat and are related each other by a scale transformation. On the 
other hand, in the case of R — > +00, our Lagrangian for the finite coupling becomes 



C 4 



1 



1 - 8e~ mR + (12mR + 28)e" 2mK (R 2 + 9 2 ) 



2mR 



(4.40) 



Consequently the equation of motion of the relative distance R for D^m) solution at the finite 
gauge coupling is found to be 



R 



-e~ mR 4m + e- 2mR 



(12m 2 R - 10m) (R 2 - 9 2 



On the other hand, the asymptotic Lagrangian for infinite coupling at R 
to be [28] 



(4.41) 

00 has been found 



£ = - m ( |i _ 4(mi? - l)e- 2rnR 



(R 2 + 9 2 ), 



(4.42) 



and the corresponding R — > 00 limit of the equation of motion of the relative distance R for 
D (m) solution at the infinite coupling reads [28] 



R 



-2mR 



(4m 2 i?-6m)(i? 2 -9 2 ). 



(4.43) 



We see that for R 2 > 9 2 there is always an attractive force operating between walls which are 
sufficiently far apart, irrespective of the strength of the gauge coupling. However, we find that 
the strength of the force for the finite coupling behaves as e~ mR , which is much larger than that 
for the infinite coupling case of e~ 2mR . This explains the reason why we obtained much stronger 
force between the walls in the case of the finite coupling. 
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4.3 Approximation in 1/g 2 and Asymptotic Behaviors 



Here we wish to discuss the previously proposed expansion in inverse powers of gauge coupling 
to obtain finite coupling results [28]. We also analyze the asymptotic behavior of the BPS wall 
solutions for generic values of gauge coupling to understand the reason why the power series 
approximation exhibits pathological behavior. 

We expand Reip in inverse powers of the gauge coupling : a = (2g 2 ()- 1 = k 2 /(Am) 2 as 

oo 

Re^ = ^VV>„- (4.44) 

n=0 

By substituting the expansion to the equation (4.5), we obtain an expansion in power of a. The 
leading power 1/a comes only from the right-hand side and gives the result at the infinite coupling 
in Eq.(4.21) : tp = W. Using this result, all the successive powers can be solved iteratively 

n=0 n=0 y n=l V ; ' \/=0 / 

Several lower order results are given explicitly by 

= W, Vi = W {2 \ i> 2 = W w + (W (2) ) 2 , 

^ = W ^ + W^W™ + 2(W^) 2 + ^(W^) 3 , 

fa = W ® + QW {%) W {2) + 12H/ (5 % (3) + 9(W^) 2 + \2W^{W {2) f 

+ 12{W^) 2 W {2) + 4(H/ (2) ) 4 , • • • , (4.46) 

where is defined by = d n W/dy n . 

First let us compare the power series approximation with our exact solution of the single wall 
62(771) for the finite gauge coupling k = 2. We illustrate the result of the power series expansion 
up to the l-th order S (/) (t/) = d y tp^(y) = J2n=o ^dy^niy) by setting a = (2g 2 (y 1 = (2m)- 2 , and 
compare it with the exact solution d y ip(y) in Fig. 9. Near the wall, the approximation oscillates 
wildly and shows no indication of convergence. 



Similarly we can obtain the power series approximations for the Kahler metric F for the 
relative position and phase for two walls. For concreteness, we will take the case of our exact 
solution D^(m) for the finite coupling k — 4, and compare the exact Kahler metric F(z)4 with 
a profile obtained by the approximation in inverse powers of the gauge coupling. We can obtain 
the approximations up to the l-th order F(z)f > , by inserting the above expansion (4.46) to the 
action (4.34) and by setting a = (2g 2 ()- 1 = m~ 2 . In Fig. 10a), we compare the power series 
approximation up to various orders starting from the zero-th order F (infinite gauge coupling 
case) up to the third order approximations F 4 . We see wild oscillations near the wall. As 
the order of approximation increases, this oscillatory behavior becomes wilder, without any 
indication of dying out. If we abandon to describe the region near the wall, we can examine 
the approximation in the asymptotic regions away from the wall. We also illustrate the Kahler 
metric F(z) for larger values of z — mR — > 00 in Fig. 10b). Even there, there is no indication 
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Figure 9: Comparison of approximations in inverse powers of gauge coupling to the exact solution 
E ex (|/) for the vector multiplet scalar E(y) = d y Reip(y) with a finite gauge coupling g 2 ( = 2m 2 . 
Approximation up to the Z-th order is denoted as E^. The zero-th order approximation (infinite 
gauge coupling) E^ is represented by a dashed line. E^ by a dash-dotted line, E^ by a short 
dashed line, and E^ 3 ) by a dotted line. The mass parameter is set to m — 1. 

for the series to converge to our exact result. Actually, we can make the discrepancy more 
quantitative. For instance let us take the power series approximation up to the fourth order 
F(mR)^ for the function F(mR) with an arbitrary coupling a = (2g 2 C,Y 1 = (k/(4m)) 2 . Its 
asymptotic behavior at R — > oo is given by 

F(mR)f = 1 + (-4(1 + Am 2 a + 16m V + 64mV + 256m 8 a 4 )mR 

+ 4 + 24m 2 a + 128m V + — mV + ^^mV + 0(a 5 )) e~ 2mR 

3 15 J 

+0(e- 4mR ), (4.47) 
whereas the exact form of the function with the finite coupling at R — > oo is 

F{mR) 4 = 1 - 8e- mR + (UmR + 28)e~ 2mR + 0{e~ AmR ). (4.48) 

It is now clear that the power series approximation F(mR)^\ I — > oo can never converge to the 
function F{mR)±. Therefore we conclude that the approximation by an expansion in inverse pow- 
ers of the gauge coupling is pathological, and may not be suitable to extract physical quantities 
correctly. 

Finally we analyze the asymptotic behavior of the BPS solution for arbitrary finite gauge 
coupling. This can be worked out even without exact solutions and may provide an alternative 
approximation scheme, rather than the expansion in inverse powers of gauge coupling. Let us 
take single wall case with two hypermultiplets for simplicity. Since BPS solutions are essentially 
real when the moduli parameters are constants (not fields), we can ignore Inr0(y) and obtain the 
BPS equation as 

d 2 ij 



= g 2 ((l-e-^{e 2my + e~ 2my )) (4.49) 
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a) 



b) 



Figure 10: Comparison of approximations in inverse powers of gauge coupling up to the Z-th 
order F±\z) to the exact solution F^z) for the Kahler metric of moduli fields as a function of 
z = mR. a) Near the region of coincident walls, b) Asymptotic region y = mR — > oo. 

Let us parametrize the gauge coupling in terms of k which is no longer an integer here : g 2 ( = 
8m 2 /k 2 . We wish to devise an approximation scheme valid for large values of z = my. Since 
the vector multiplet scalar X = d y ip(y) should approach to a constant asymptotically, the BPS 
equation dictates that i]){y) should approach ±my at y — > ±00. We expand the remaining 
subleading order terms in a series which will be determined successively as solutions of iterative 
equations 



where we can assume that the remaining terms should vanish asymptotically, since we can absorb 
a possible constant term into a shift of y 




(4.50) 



n=l 



e n (y)^0 asy^O, e^y) > e 2 (y) > • • • . 
By substituting the expansion to the BPS equation, we obtain 



(4.51) 




(4.52) 



We can rewrite this equation with z 



my and g 2 ( = 8m 2 /k 2 




n=l 




(4.53) 



The first order approximation is determined by the following equation 




(4.54) 



22 



which yields a general solution with an arbitrary integration constant c\ 



4z 



1 




(4.55) 



e 1 = 2cie k + 



where we have used the boundary condition E\ — > at y — > oo to eliminate another possible 
integration constant. Therefore the asymptotic behavior of the first order approximation e\(y) 
depends on the strength of the gauge coupling constant g 2 = j£ as 



The wall solution requires that the vacuum value m is reached from below : E(y) = d y tp(y) < m 
as y — > oo, This implies £\{y) > 0. Therefore we obtain up to the first order approximation 



This result can be regarded as the origin of a different asymptotic behavior for the finite coupling 
exact solution (k — 2 case) compared to the infinite coupling solution (k — case). It is 
interesting to observe that the critical case of k — 1 has an exceptional asymptotic behavior 
corresponding to degenerate exponents. This complication may be related to the fact that we 
have not yet succeeded to obtain an exact solution in that case. 

It is now straightforward to extend our approximation scheme to more general multi-wall 
cases. Since our approximation scheme from asymptotic region is applicable to any values of 
gauge coupling constant, this should be useful as an alternative approximation method instead 
of the expansion in inverse powers of gauge coupling. Our iterative approximation is at least 
guaranteed to be valid for asymptotic region y — > oo, although it may wildly oscillate and possibly 
diverge near the wall. On the other hand, we have seen that the expansion in power of inverse 
gauge coupling does not converge near the wall nor asymptotically. To obtain a full solution, one 
needs to determine the integration constant such as C\ which can be done by smoothly connecting 
to solutions from the region near the wall. Only when this integration constant is chosen to be 
a particular value, the solution should smoothly connect to the solution near the wall. 
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k > 1, ci > 
k = 1 
< k < 1 



(4.56) 




k > 1 
k= 1 
< k < 1 
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Abstract 

A series of exact BPS solutions are found for single and double domain 
walls in M = 2 supersymmetric (SUSY) QED for finite gauge coupling 
constants. Vector fields are found to be massive, although it is localized 
on the wall. Massless modes can be assembled into a chiral scalar mul- 
tiplet of the preserved M = 1 SUSY, after an appropriate gauge choice. 
The low-energy effective Lagrangian for the massless fields is obtained 
for the finite gauge coupling. The inter-wall force is found to be much 
stronger than the known infinite coupling case. The previously proposed 
expansion in inverse powers of the gauge coupling has pathological oscil- 
lations, and does not converge to the correct finite coupling result. 
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1 Introduction 



In recent years, much efforts have been devoted to the brane world scenario [?]-[?], where our 
world is supposed to be realized on a subspace such as a wall embedded in higher dimensional 
spacetime. On the other hand, supersymmetry (SUSY) has been one of the most fruitful ideas to 
build realistic unified models beyond the standard model [?] . SUSY helps to obtain domain walls 
and other stable solutions by requiring a part of SUSY to be preserved by the field configuration. 
Then the configuration has minimum energy with the given boundary condition and is automat- 
ically a solution of field equations. Wall configurations can preserve half of SUSY charges [?]-[?], 
which are called | BPS states [?]. To obtain a wall with the four-dimensional world volume, we 
need to consider a fundamental theory in five or more spacetime dimensions. If such a theory are 
supersymmetric, it must have at least eight SUSY charges [?]. Simplest theories with eight SUSY 
consist of hypermultiplets. The nontrivial interactions of hypermultiplets require either nonlin- 
earity of kinetic term (nonlinear sigma model) or gauge interactions [?]. The target spaces of 
such nonlinear sigma models are hyper-Kahler (HK) manifolds [?, ?] and can most conveniently 
be obtained by using vector multiplets as Lagrange multiplier fields [?, ?]. By introducing the 
Fayet-Iliopoulos (FI) term for the U(l) Lagrange multiplier vector multiplet and mass terms for 
hypermultiplets, one can obtain nontrivial potential terms with N discrete SUSY vacua for N 
hypermultiplets. Domain wall solutions have been obtained in these massive nonlinear sigma 
models with the target space T^Qp-W- 1 [?]-[?] and their generalizations [?]. One of the most 
convenient methods to obtain these nonlinear sigma models is to take a limit of infinite gauge 
coupling from the usual minimal gauge interactions of hypermultiplets, namely massive M = 2 
SUSY QED. However, it has been known that the massive Af = 2 SUSY QED itself has discrete 
SUSY vacua and may allow domain walls even with finite gauge coupling [?]-[?]. 

The massive nonlinear sigma models with multi-flavor have been found to give BPS multiple 
wall solutions containing a number of parameters, which represent the relative positions and 
phases of these multiple walls and are called moduli [?]. When these moduli parameters are 
promoted to functions of the world volume coordinates of the walls, one can obtain the low- 
energy effective Lagrangian of moduli fields, which are massless [?, ?, ?]. It has been noted 
that gauge field should acquire nontrivial configuration when these moduli depend on the world 
volume coordinates [?]. It has also been proposed that an expansion in inverse powers of gauge 
coupling may be useful to obtain the moduli dynamics at finite gauge coupling [?]. 

To fulfill the goals of the brane-world scenario, we need to localize particles of the standard 
model to "branes" , such as domain walls. It has been a notoriously difficult problem to obtain 
massless localized gauge field on a wall [?, ?, ?]. The general arguments [?, ?] and an explicit 
model [?] suggest that gauge fields may be localized on a wall, but will not be massless in the 
massive Af = 2 SUSY QED with finite gauge coupling. Nevertheless, it is still interesting to 
study explicitly the role and the properties of dynamical gauge fields in the interacting system 
of the vector multiplet and hypermultiplets. In this respect, it should be useful to obtain exact 
solutions with finite gauge couplings [?] 1 . In the four-dimensional massive Af = 2 SUSY QED 
which can be obtained by a dimensional reduction of our model, it has been argued that a massless 
gauge field is obtained by dualizing the massless Nambu-Goldstone scalar [?] 2 . While this is an 

1 In Ref.[?], an exact solution of domain wall junction has been obtained for the massive M — 2 SUSY QED 
with the finite gauge coupling for the first time. Therefore it automatically implies an exact solution of domain 
wall for finite gauge coupling. 

2 They are interested in an Af = 1 gauge theory obtained from Af = 2 SU{2) gauge theory deformed by an 
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intriguing result, the gauge field as a dual of a compact scalar is an intrinsically three-dimensional 
peculiarity which is not valid to the more realistic situation such as our five-dimensional theory. 
In another paper, we show that massless localized gauge fields can be obtained by coupling our 
model to a tensor multiplet [?] . 

The purpose of our paper is to work out exact solutions of BPS single and multiple walls 
with finite gauge coupling in massive M = 2 SUSY QED in five dimensions and to study the 
fluctuations on the domain wall background. We obtain a series of exact solutions for single as 
well as multiple walls provided gauge coupling has particular values relative to the ratio of the 
mass parameters and FI parameters. We also show that the gauge multiplet is massive for any 
finite gauge coupling, irrespective of the possible existence of exact solutions. This is consistent 
with the previous general arguments [?, ?] and explicit examples[?]. In our model, the massless 
Nambu-Goldstone scalar still exist, but are apparently unrelated to the dual gauge field. Instead, 
we obtain the Nambu-Goldstone scalars and the Nambu-Goldstone fermions which form a chiral 
scalar multiplet of the remaining Af = 1 SUSY, by choosing an appropriate gauge. Moreover 
we have only massive spectra for the vector field. These results are in contrast to the previous 
results of the massive M = 2 SUSY U(l) gauge theory in four dimensions, where the Nambu- 
Goldstone scalar is identified as the dual of the three-dimensional gauge field [?]. As a result 
of exact solutions, we also obtain the low-energy effective Lagrangian for massless moduli fields 
for finite gauge couplings, by promoting the moduli parameters into four-dimensional massless 
fields [?, ?, ?]. In the process of carrying out the program, we observe that the gauge field should 
acquire nontrivial configuration with nonvanishing field strength around the wall, confirming the 
observation in Ref. [?]. Moreover, we find that this nontrivial gauge field configuration should 
be determined by means of the field equation of gauge field. We compute also some lower order 
terms of the previously proposed expansion in inverse powers of gauge coupling [?], and compare 
the expansion with our exact result. We find that the expansion oscillates wildly and does not 
converge to the exact result in any smooth or uniform way. Therefore the expansion has a 
difficulty or at least subtlety in extracting physical quantities. 

In Sect. 2, we introduce our model of massive M = 2 SUSY QED and obtain its SUSY vacua 
and BPS equations. In Sect. 3, we study the single wall in the case of two hypermultiplets 
(N = 2), and work out fluctuation spectra to show that there are massless Nambu-Goldstone 
modes, but that the vector field has only massive fluctuations. In Sect. 4, a series of exact 
solutions are constructed for N — 1 walls in the case of N hypermultiplets and obtain the low- 
energy effective Lagrangian for the massless fields. We also work out a few lower order terms of 
the expansion in inverse powers of gauge coupling and compare it to our exact solution. 

2 Massive M = 2 SUSY QED and BPS equations 

The simplest building block of M = 2 SUSY theory in five dimensions is hypermultiplets which 
consist of SU(2)r doublets of complex scalar fields H tA , Dirac fields ip A and complex auxiliary 
field F A , where i — 1,2 stands for SU (2) R doublet indices and A = 1, - ■ • ,N stands for flavours. 
For simplicity, we assume that these N hypermultiplets have the same U(l) charge, say, unit 
charge. The U(l) vector multiplet consists of a gauge field W M) a real scalar field S, SU(2) R 



adjoint scalar mass term. However, they observed the model reduces to our Af = 2 SUSY QED to the leading 
order of the mass deformation parameter. 
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doublet of gauginos A* and SU (2) R triplet of real auxiliary fields Y a , where M, N = 0, 1, • • • , 4 
denote space-time indices, and a = 1,2,3 denotes SU(2) R triplet index, respectively. In this 
work, we shall consider a model with the minimal kinetic term for hypermultiplets and vector 
multiplets. The M = 2 SUSY allows only a few parameters in our model : the gauge coupling 
g, the mass of the A-th hypermultiplet m A , and the FI parameters ( a . The FI parameters are 
real and transforms as a triplet under SU(2) R . Then the bosonic part of our Lagrangian reads 

[?]"[?],[?] 

£boso„ = -^- 2 {F MN {W)f + -^(9 M S) 2 + (V M H)\ A (V M H iA ) - Hl A {H - m A fH iA 
^9 ^9 

+7^(Y a ) 2 - CY a + H\ A {a a Y a ) i jH jA + F A F A : (2.1) 

29 

where a sum over repeated indices is understood, Fmn(W) = OmWn — 9nWm, covariant deriva- 
tive is defined as Dm = 9m + iWm, and our metric is i]mn = (+1, — 1, • • • , — 1)- 

To ensure the discreteness of SUSY vacua, we need to make all hypermultiplet masses non- 
degenerate. Without loss of generality, we assume the following ordering of the hypermultiplet 
mass parameters 

m A+1 < m A (2.2) 

for all A. When all hypermultiplet masses are nondegenerate, the symmetry of our model reduces 
to U(1) N . The diagonal U(l) is gauged by the vector multiplet W M) and the remaining U(l) N ~ l 
is the global symmetry. 

The easiest way to find SUSY vacua is to explore the condition of vanishing vacuum energy. 
To facilitate the procedure, let us first write down equations of motion of auxiliary fields Y a and 
Ft 

Y a = g 2 [C-H\ A (o a y 3 Hi% (2.3) 
F A = 0. (2.4) 

After eliminating the auxiliary fields, we obtain the potential V pot 

£boso„ = ~(F MN (W)) 2 + ^(d M E) 2 + {V M H)\ A {V M H iA ) - V^, (2.5) 
^9 ^9 



v - g2 



C-H\ A {a a y 3 W A \ +H\ A {Y,-m A ) 2 H iA . (2.6) 

The vanishing vacuum energy is achieved by requiring 

HU^YjH^ = C, (2.7) 
(E - m A )H iA = 0, (2.8) 

where, in the second equation, index A is not sumed. These SUSY vacuum conditions guarantee 
the full preservation of SUSY and can also be derived by requiring the SUSY transformation of 
fermions to vanish as we see immediately. 

By making an SU(2)r transformation, we can always bring the FI parameters to the third 
direction without loss of generality 

C = (0,0,C), C>0. (2.9) 
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In this choice, we find N discrete SUSY vacua (A = 1, • • • , N) explicitly by solving (2.7) as 

S = m A , \H 1A \ 2 = (, H 2A = 0, 

H 1B = 0, H 2B = 0, (B^ A). (2.10) 

Since fermions are assumed to vanish in the wall configuration, we need to examine only 
SUSY transformations of fermions to find a configuration preserving a part of SUSY. Gaugino 
A* and hyperino ip A transforms as 3 

5 £ \ l = (^7 M ^M7v(lU)+7 M 9Ms)£ l + ^YV a )V, (2.11) 

5 e ^ A = -iV2~U M V M H lA + i(E - m A )H iA \ e^eP + V2F A e\ (2.12) 

To obtain a wall solution, we assume the configuration to depend only on the coordinate of one 
extra dimension, which we denote as j/ = x 4 . We also assume the four-dimensional Lorentz 
invariance in the world volume coordinates ), which implies 

F MN (W) = 0. (2.13) 

Since we are interested in the \ BPS configuration, we require the above SUSY transformations 
(2.11) to vanish for half of the Grassman parameters specified by 

P+e 1 = 0, P_e 2 = 0, (2.14) 

where P± = (1 ± 7s)/2 are the chiral projection operators. Finally we need to eliminate the 
auxiliary fields Y a and F A by their algebraic equations of motion (2.4) to make the BPS condition 
as first order differential equations for physical fields. Thus we obtain the \ BPS equations for 
the massive M = 2 SUSY QED as 

dyV = g 2 (C - H\ A H 1A + H\ A H 2A ) , (2.15) 

2g 2 Hl A H 1A = 2g 2 H\ A H 2A = 0, (2.16) 
V y H iA = (m A -Z)H iA , i = l,2, A=1,---,N. (2.17) 

One can easily see that the translation invariant vacuum requires the vanishing of the left-hand 
side of Eqs.(2.14) and (2.16), which implies the same condition as the full preservation of SUSY 
in Eq.(2.7). 

We are interested in the solution of these BPS equations which interpolate two different vacua 
in Eq.(2.9). Since the BPS equation is a first order differential equation, the boundary condition 
dictates that H 2A should vanish identically 

H 2A (y) = 0. (2.18) 

The energy density of the BPS solution can be found by making a Bogomolny completion of 
the energy functional 

T w = dy{^- 2 (d y J: - g 2 (( - \H 1A \ 2 )) 2 +\V y H 1A - (m A - S)iJ^| 2 } 



+ 



1 y=+oo 

(Z + (m A -Z)\H 1A \ 2 . (2.19) 

; -I i/=-00 



3 Our gamma matrices are 4x4 matrices and are denned as : {7 M , 7^} = 2r\ MN , ^ MN = ^ [y M ,^ N ] — j\- m j n ^ , 
y 5 = Z7°7 1 7 2 7 3 = — ij 4 . 
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yielding the BPS wall tension for solutions interpolating the A-th vacuum and the _B-th vacuum 

T w = ((m A - m B ), (2.20) 

assuming A > B. The structure of the above BPS equations shows that only those hypermultiplet 
scalars H 1C ^ with A > C > B have nonvanishing values besides the vector multiplet scalar 
S [?], [?]. Defining I = A — B, we call such a BPS configuration as / wall solutions, since it 
represents / separate walls at least when these walls are sufficiently far apart, as we shall see in 
the Sect. 4. 



3 Single wall BPS solutions 

3.1 Exact solutions for single wall 

Since single wall interpolates two adjacent SUSY vacua, we shall assume N = 2 (two hypermulti- 
plets), without loss of generality. For simplicity, we take mass parameters of the hypermultiplets 

as 

mi = — m 2 = m. (3.1) 
The boundary conditions for BPS equations are given by 

£(— oo) = —m, S(oo) = m, (3-2) 
tf n (-oc)=0, tf n (oo) = v/C, 

H l2 (-oo) = y/C, H 12 (oo) = 0. (3.3) 

The set of the BPS equations (2.14)-(2.16) with the above boundary conditions (3.2)-(3.3) 
are known to be solved exactly for infinite gauge coupling [?], [?]. Recently another exact solution 
has been found for finite gauge coupling, provided the gauge coupling g satisfies the following 
relation with the ratio of the FI parameter ( and the mass difference squared (2m) 2 [?] 

g 2 ( = 2m 2 . (3.4) 

By generalizing this exact solution for finite gauge coupling, we find that there are exact solutions 
for the following values of gauge coupling g 

i2 



9X = -p, (3.5) 

with appropriate integers k. We shall denote the single wall solution for k as Sk{m). The integer 
k indicates the value of the gauge coupling g through the relation (3.5), and the mass parameter 
m indicates the tension T w of the single wall in unit of the FI parameter as T w = 2(m. The 
infinite gauge coupling corresponds to the case k — 0, and the previously obtained finite gauge 
coupling case to k — 2 [?]. We have explicitly obtained exact solutions in the case of k — 0, 2, 3, 4. 

In this section, we will describe the exact solution S 2 (m) for k = 2 as the simplest example, 
which plays an important role in our paper on massless localized vector field [?] : 

= mtanhmy, (3-6) 

H 11 (y) = -^e +my sech my, (3.7) 



V2g 
m 

VTg ( 



771 

H l2 (y) = ^e- my sechmy. (3.8) 
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Let us recall that we can introduce two arbitrary integration constants in these solutions. We can 
have y — > y — yo corresponding to the spontaneously broken translation invariance, and we can 
also have the multiplication by a phase e~ ia ( e ia ) for hypermultiplets if 11 (if 12 ) corresponding 
to the spontaneously broken global U(l) invariance. They constitute two moduli of the solution. 
The yo has a physical meaning of the position of the wall. The vector multiplet scalar £ and 
the hypermultiplet scalars if 11 and H 12 of the solution S^m) are illustrated as a function of the 
coordinate y in the extra dimension in Fig.l. We will present other exact solutions of single wall 
after introducing a slightly different notation suitable also for multiple walls. 



S/m V2gH 1A /m 




a) b) 

Figure 1: The exact BPS wall solution S2 (m) for finite gauge coupling g = \f2mj\fC, and the 
tension T w = 2m(. a) Scalar field E(?/) of vector multiplet in Eq.(3.6) divided by the mass 
parameter m as a function of the coordinate y times m. b) Hypermultiplet scalar field H ll (y) 
in Eq.(3.7) as a function of the coordinate y times m (solid line), and H l2 {y) in Eq.(3.8) as a 
function of the coordinate y times m (dotted line). 



3.2 Fluctuation around the BPS single wall background 

Let us now turn our attention to spectra of fluctuations around the BPS single wall background. 
These fluctuation fields can be expanded into modes defined on the background. Among various 
modes, we are particularly interested in two kinds of modes : massless modes (zero modes) 
and fluctuation of vector fields. The former modes play essential role in discussing dynamics 
on the background at low-energies. The latter is to examine the role played by the dynamical 
vector field at finite coupling, rather than the role as the Lagrange multiplier field in the case of 
corresponding nonlinear sigma models. 

Zero modes usually arise from several reasons. The first one is the Nambu-Goldstone modes 
as a result of the spontaneous breaking of continuous global symmetry. The second possibility is 
the result of remaining symmetry, such as the conserved SUSY [?]. The third possibility is the 
result of index theorems, such as fermion zero modes around instantons and other solitons [?]. 

To illustrate the method, we shall work on the fluctuations on the background of the simplest 
exact solution S^m) in Eqs.(3.6)-(3.8). In the case of a single wall, such as S^m), the only 
spontaneously broken global internal symmetry is the U(l) symmetry which rotates the phase 
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of hypermultiplets oppositely 

There are two bosonic global symmetry which are spontaneously broken : translation and the 
global U(l) in Eq.(3.9). The mode function of the Nambu-Goldstone boson corresponding to the 
spontaneously broken translation is given by differentiating the background field configuration 
with respect to the coordinate in extra dimension y. We shall denote the corresponding four- 
dimensional effective field as Re0 o (^)- Similarly the Nambu-Goldstone boson associated with the 
global U(l) symmetry is given by acting an infinitesimal U(l) transformation. We denote the 
corresponding four- dimensional effective field as Im0o(^)- By writing out only zero modes, we 
thus obtain scalar fields of vector multiplet E(x, y) and that of hypermultiplets H 1A (x, y),A= 1,2 

E(x,y) = (£(y)> - d y {E(y))Re<f> (x), (3-10) 
H n (x,y) = (H n (y)) - d y (H n (y))ReMx) ~ *MH U (y))ImMx), (3.11) 
H 12 (x,y) = (H 12 (y))-d y (H 12 (y))ReMx)+im(H 12 (y))lmMx), (3.12) 

where we multiplied the factor m in front of the £7(1) Nambu-Goldstone mode wave function for 
later convenience, and the (• • •) denotes the background field configuration. These configurations 
together with vanishing gauge field Wm = are solutions of the BPS equation if Re0o and Im0o 
are constants (moduli) and hence automatically satisfy the field equations. However, we are 
now interested in the case where these Re0 o and Im0 o are four- dimensional fields, which depend 
on the coordinate x M of the world volume. Then these field configurations are no longer BPS 
configurations, nor satisfy the field equations. Since four-dimensional effective fields are defined 
as solutions of the field equations, we should demand these field configurations to satisfy the 
linearized field equations. Since Wm must vanish for constant <po, it should be proportional to 
d^<f)o(x). Therefore the Lorentz invariance in the four-dimensional world volume dictates 

W,(x,y) = 0, (3.13) 

and 

W lt (x,y)=w(y)d l Jm<f>o(x). (3.14) 

with some function w (y) to the linearized order which we are interested in. Then the linearized 
field equations for the vector multiplet scalar S, hypermultiplet scalars H lA are satisfied by 

d^Mx) = 0, (3.15) 

which just means that these Nambu-Goldstone fields are massless. 

The remaining linearized field equation for the gauge field Wm reads 

d v d v W, - d 2 W, - d,{d v W v ) + d,d y W 4 = -2g 2 ((H 1A )d,lmh 1A + \{H 1A )\ 2 W,) , (3.16) 

dvVWt - d y {d v W v ) = 2g 2 (d y (H 1A )lmh 1A - (H 1A )d y lmh 1A - \(H 1A )\ 2 W±) . (3.17) 

The solutions for Wm in Eqs.(3.13) and (3.14) are enough to satisfy the field equation for W4 in 
Eq.(3.17). The remaining field equation (3.16) for W^ determines the y dependence w(y) of the 
gauge field W^(x,y) in Eq.(3.14) as 

Wn(x,y) — (T,(y))dfjlm(f) (x). (3.18) 
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This solution shows that the corresponding field strength is nonvanishing around the wall 

F^{W){x,y) = —d y {E(y)) c^Im^o (x) . (3.19) 

In a similar model in four-dimensions, it has been observed that the gauge field must be nontrivial 
and field strength must be nonzero around the wall if the moduli depends on world volume 
coordinates [?, ?]. Our result in five dimensions are in agreement with their observation. Our 
new point is perhaps that the necessary gauge field configuration can be determined explicitly 
by using the linearized field equations for gauge fields. 

SUSY is another global symmetry which is spontaneously broken. We can obtain the Nambu- 
Goldstone fermion corresponding to the broken SUSY by evaluating the SUSY transformations 
of fermions on the background. One set of ££/(2)-Majorana spinor X l ,i = 1,2 in five dimensions 
can be decomposed into two Majorana spinors A+, A_ in four dimensions as 

A + = P+A 1 + P_A 2 , A_ = -P_A X + P+A 2 . (3.20) 

The SUSY transformations in Eqs.(2.10) and (2.11) can be evaluated by using the BPS equations 
(2.14)-(2.16) to yield for the gauginos 

5 £ X + = i [d y (E) +g 2 {C~ \(H 1A )\ 2 )} (P + e + - Ps + ) = 2id y {T) (P + e + - P_s + ) 

5 £ \- = i [d y {T) -g 2 {(- \(H 1A )\ 2 )} (P_e_ + P + e_) = 0. (3.21) 

Since e_ is preserved, it drops out from these transformations. If we evaluate explicitly for our 
exact solution £2(771) in Eqs.(3.6)-(3.8), we obtain 



5 eK = .._ l2 — (P+£+ ~ P-e+) , S e X- = 0. (3.22) 



2im 2 
cosh 2 my 

Similarly the hyperino component of the Nambu-Goldstone fermion wave function is given by 

6 e if> A = V2 [d y (H 1A ) + ((E) - m A ) (H 1A )] P + S- + V2 [-d y (H 1A ) + ((E) - m A ) (H 1A )] P_e + 
= -2V2d y (H 1A }P_E + . (3.23) 

For the exact solution S 2 (m) in Eqs.(3.6)-(3.8), it becomes explicitly as 

5 e ^ = -5 £ ^ 2 = -——^P^e + . (3.24) 
g cosh my 

The above transformation laws (3.21) and (3.23) show that — 4e + is proportional to the zero 
momentum component of the Nambu-Goldstone fermion. Therefore we can define the Nambu- 
Goldstone fermion field Xo( x ) °f the four- dimensional effective low-energy theory on the world 
volume of the wall as 

\ + (x,y) = - l -d y (X(y)) l5X o(x), \-(x,y) = 0, (3.25) 
^ A (x,y) = ^=d y (H 1A (y))P^ X o(x), (3.26) 

suppressing to write massive modes. We see explicitly that the Nambu-Goldstone fermion appears 
only with the left-handed chirality in the hypermultiplets as is usually dictated by index theorems 
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for fermions localized on domain walls [?]. We can also verify that the Nambu-Goldstone fermion 
satisfies the linearized equations of motion with a vanishing mass eigenvalue, by using the BPS 
equations (2.14)-(2.16). 

Now let us consider the requirement of the preserved symmetry, especially SUSY. The SUSY 
transformation property under the preserved SUSY specified by Eq.(2.13) is given by 



S e _W lt = -ie^X-, (3.27) 

S £ _W 4 = e_A+, (3.28) 

5 £ _S = -i£_7 5 A+, (3.29) 

5 £ _H 1A = -V2e-P-4> A , (3.30) 

5 £ _H 2A = -V2e_P + ^ A . (3.31) 



This transformation property does not fit well with the above massless particles of Nambu- 
Goldstone modes in Eqs(3.10)-(3.14), (3.25), and (3.26), especially for the vector multiplet. 
However, we should remember that we have a freedom to make gauge transformations to make 
the above massless particles in conformity with the SUSY transformation properties. Let us 
perform the following gauge transformation, which is proportional to the fluctuation field Im0 o 

W 4 — > W 4 -9„«S)Im0o), 
W, — W M -^((S)Im0 o ), 
H 1A — > (H 1A ) + i ((S)lm0 o ) (H 1A ). (3.32) 

Then the hypermultiplet scalars become 

H u (x,y) = (H u (y)) - ^(/J 11 (y))Re0 o (x) -i(m- (E(y)» (H u (y)) lmfo(x) 

= (H u (y)) - d y (H u (y)) (Re0 o (:r) + Tm0 o (x)) , (3.33) 

H 12 (x,y) = (H 12 (y)) - ^(iJ 12 (y))Re0 o (x) -i(-m- (E(y)» (H^^lmfa (x) 

= (H 12 (y)) - d y (H 12 (y)) (Re0 o (^) + Tm0 o (x)) . (3.34) 

Therefore we find that the two Nambu-Goldstone bosons corresponding to translation and £7(1) 
global rotation forms a complex scalar (f>o(x) of a chiral scalar multiplet of the preserved four 
SUSY e_. At the same time, the four-dimensional gauge fields vanishes, and vector multiplet 
scalar E and the fourth (extra dimension) component of vector field W± fit into the complex 
scalar <po(x) 

Z(x,y)+iW 4 (x,y) = (E(y))-^(E(y))(Re0 o ( a ;)+Tm0 o ( a ;)), (3.35) 
W li (x,y) = 0. (3.36) 

These mode expansions together with the fermionic ones in Eqs.(3.25) and (3.26) are now consis- 
tent with the preserved SUSY which defines the four-dimensional M = 1 SUSY transformations 
for the massless fields in the effective field theory 

8 e _Mx)=e-P-Xo(x)- (3.37) 

Thus we find that all the Nambu-Goldstone bosons and fermions together form a chiral scalar 
multiplet (<f>o(x) , Xo( x )) °f ^ ne preserved SUSY. Moreover we do not find any index theorem to 
force us additional zero modes. Therefore we do not expect any more zero modes. 
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Let us now turn to the issue of mass spectra, especially for vector fields. The linearized 
equation of motion for vector field and W 4 are given in Eqs.(3.16) and (3.17). By exploiting 
the U(l) gauge invariance, we can always choose a gauge where W±{x,y) = identically. We can 
also decompose vector field into the transverse part and longitudinal part w 

W^x : y) = W fl (x,y) + d ll w(x,y), = 0. (3.38) 

The transverse part of the vector field 4 is decoupled from the rest of the fluctuations because 
of the four-dimensional Lorentz invarince and gives the following mode decomposison with the 
mode function a n (y) for the mass eigenvalue in n of the four-dimensional effective field wjP\x) 

W ll (x,y)=Y i a n (y)Wlr\x), (3.39) 

n 

-d 2 a n (y) + V(y)a n (y) = m 2 n a n (y). (3.40) 

where the potential V(y) is defined by 

V(y) = 2g 2 \(H 1A )\ 2 : (3.41) 

and is illustrated in Fig.2 for the background S 2 (m) in Eqs.(3.6)- (3.8), and S 3 (m) and S^m), 
whose exact solutions are given explicitly in Sect. 4.1. These potentials have a dip near the 




Figure 2: The potential V(x) of Eq.(3.41) for ^(m) (solid line), S^m) (dotted line), and S 4(171) 
(dashed line) in the eigenvalue equation for the vector fluctuation (3.40) as a function of extra 
dimension coordinate y. The mass parameter is taken to be m = 1. 



wall and approaches 2g 2 ( at y — > ±00. Therefore all the vector fluctuations become infinitely 
heavy as g — > 00. Potentials are always positive definite, but have attractive forces due to a dip 
around the wall for any finite gauge coupling. It is likely that there may be discrete low mass 
states before reaching a continuum starting from m 2 = 2g 2 (. The potential becomes infinite for 
the infinite gauge coupling. It is clear that the eigenvalue can never vanish : m n > 0. It is also 
easy to see that m n are of the order of m for, say, S^to). This result is consistent with the 

4 This component corresponds to the polazation states proportional to the polazation vector e l ^,i — 1,2,3 
defined by e^p^ = in momentum space, representing a genuine four-dimensional vector field. 
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general argument that the gauge field will obtain a mass of the order of the inverse width of the 
wall. Although the low-lying vector modes have presumably discrete spectra and are localized 
around the wall, all these modes are massive. To obtain a massless vector field, we can add 
tensor multiplet [?]. 

To find out the remaining field equations, we define fluctuations of field as E = (E) + s, 
H Al = (H At ) + h Al . The linearized field equations are given by 

= d M d M s + 2g 2 {\H Ai \ 2 )s + 2<7 2 /4<(E - m A )H M ) + 2« 7 2 (4 4 (E - m A ))h A \ (3.42) 

= d M d M h Al + i(d M W M )(H Ai ) + 2iW M d M (H Ai ) - g 2 (( - (\H Bj \ 2 ))h Ai 

+ g 2 (H Ai ) ({H^)h B i + (/F B )/4) + 2((E) - m A )(H Ai )s + ((E) - m A ) 2 h A \ (3.43) 

These linearized field equations can be decomposed into three sets of coupled equations. The 
first set consists of the fluctuations of vector multiplet scalar E, and those of the real part 
of the hypermultiplet scalars KeH , A = 1,2. The second set consists of imaginary part 
of the hypermultiplet scalars lmH 1A , A = 1, 2 and the longitudinal part of the vector field 
w in Eq.(3.38). The third set consists of the lower components of the hypermultiplet scalars 
H 2A , A = 1,2. Among various massive modes, we find the following tower of massive modes as 
a solution of the first set of coupled equations 

s = Y,( d v a n{y))<t>n{x), h A1 = - ^(iJ A1 (t/))a n (t/)0„(x), (3.44) 

n n 

with h A2 = Wm = 0. This effective field 0„(x) has the identical mass squared m 2 as the effective 
field Wj?\x) of the transverse vector in Eq.(3.40). Moreover the mode function is also given by 
the same wave function a n (y) as the transverse vector. Therefore it is likely that these effective 
fields are related by the preserved symmetry, in particular the Af = 1 SUSY. In fact, a massive 
vector multiplet of J\f = 1 SUSY should contain a scalar particle. 



4 BPS solution for Multiple walls 

4.1 Exact solutions with finite gauge couplings 

It is convenient to introduce a complex function ip(y) to solve the BPS equation for hypermulti- 
plets (2.16) [?] 

N-2 

H 1A {y) = ^/(exp(-i>(y)+m A (y-y ) + J2aA a ra), (4.1) 

a=l 

where the complex parameters r a , a — 1, • • • , N—2 are collective coordinates arising as integration 
constants. Since two complex parameters among the integration constants can be absorbed by 
a shift of tp and y , we can choose the N x (N — 2) fixed real matrix a A a to be of rank N — 2. 
The BPS equation for hypermultiplets is equivalent to the following equation for ip(y) 

d y ^ = Z + iW 4 . (4.2) 
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Because of the vanishing field strength (2.12), the vector field W 4 = lmd y ip is a pure gauge. 
However, when we consider dynamics of domain walls, this term will play an important role. 
We have the following moduli parameters in the solution (4.2), 

yo = Y + i9 , (4.3) 
r a = R a + i6 a , (4.4) 

with one and N — 2 complex dimensions, respectively. The Y and R a are related to the center 
of mass and relative positions of the a-th domain wall, respectively, and 6q and 8 a to the overall 
phase and relative phases of the a-th wall, respectively. Using the variable ip, the BPS equation 
for vector multiplet (2.14) becomes 

j N N-2 

— d 2 y Re(4>) = 1 - ^exp(-2Re(V) + 2m A {y - Y ) + 2 ^ a A a R^j 

^9 A=l a=l 

= 1 - exp (-2Re(V>) + 2W) . (4.5) 
Here, the explicit dependence on y in Eq.(4.5) can be assembled by defining a function W(y) as 

N N-2 

W = log exp (2m A (y - Y ) + 2 ^ au^a) . (4.6) 

A=l a=l 

To ensure that d y T, = at y — > ±oo, we should impose boundary conditions for the above 
equation as 

AT-2 



Re(^) mi( 2 /-y ) + 53ai°i2 a , y -> +oo, (4.7) 

o=l 
JV-2 

Re(^) — > - F ) + ^ ajv^a, y -> -oo. (4.8) 



a=l 



To define an appropriate variable for the position of each wall, let us consider a configuration 
where only two adjacent exponential terms in the sum of Eq.(4.5) are large and the others are 
negligible, then the function W has a profile for a single wall. Therefore it is natural to define 
the position y A of the A-th wall as 



(AT-2 \ / N-2 

2m A (y A - Y ) + 2 Y ^A a R a = exp 2m A+1 (y A - Y ) + 2 ^ a A+1 a R a 
a=l J \ a=l 

-+ y A = Y - ^^ a ' aAa)Ra . (4.9) 

m A+1 - m A 

This moduli parametrization y A has an intuitive meaning of the position of the A-th wall and 
Ua — Ua+i corresponds to a distance between two walls, at least when the distance is large. 
Relations between y A and the relative positions R a defined in Eq.(4.4) are obtained by 

N-2 A-l 

- J]K a - «i a R = Yl ( m ^+i - (Vb -Y ), A = 2, • • • , N. (4.10) 

a=l B=l 
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By choosing A = N, we find 

~N-1 N-2 

^2 ( m B+i - m s) Vb + ^2(a a N - al)R a 



rriN — mi 



B=l o=l 



(4.11) 



which becomes the center of mass coordinate when a a N = a". Because of the translational 
invariance, Re(^) should be a function of real variables y — Y and R a . From now on, we will 
take Y = unless otherwise stated. 

Now we will present a series of exact solutions for single and double walls for finite values 
of gauge coupling. To compare single and multi-wall solutions, let us assign the following mass 
parameters for hypermultiplets for double walls (N — 3) : 

m A = {m,0, -m). (4.12) 

This convention is intended to make the total energy density of the double wall to be identical 
to the single wall, so that the double wall situation can be most naturally compared to the single 
wall situation with the same energy density (tension) in the coincident limit of two walls, since 
the mass parameters of the single wall is assigned to be m A = (m, —m) in Eq.(3.1), and the total 
energy density is just given by the difference of the two extreme masses as given in Eq.(2.19). 
For the moduli parameter, we use y\ — y 2 = R 1 = R and choose the a a 1 as 

a/ = (0,m/2,0), 2a A a R a = 2a A 1 R 1 = (0, mR, 0), (4.13) 

since the rank of the matrix aA a is N — 2. This relative distance appears only in multiple wall, 
but not in the single wall. 

The function W = W S i Qg i e for the single wall case becomes 

W single ^log(e 2 ^ + e- 2ms 0. (4-14) 
Similarly the function for the double wall case W = Wdouhie is found to be 

WWbie = \ log(e 2m f + e~ 2m y + e mR ). (4.15) 

The solvable cases of finite gauge coupling are found to be 

8m 2 

T 2 

where the mass parameter m is defined for single wall case in Eq.(3.1), and double wall case 
in Eq.(4.12). We will denote the single wall solution as Sk{m) and the double wall solution as 
Dk(m) with the coupling defined by (4.16) with k and the mass parameter m. 

Let us list the exact solutions that we are able to obtain 



9 2 ^—, k = 0,2, 3, 4, (4.16) 



SoM 


Reifj = 


hog{e 2my + e- 2my ), 


(4.17) 


S 2 (m) 


Reip = 


log(e my + e- my ), 


(4.18) 


S 3 (m) 


Re^ = 


^log(e§^ + e-§^), 


(4.19) 


S A {m) 


Re^ = 


log(e my + e~ my + 


(4.20) 


D (m) 


: Re^ = 


Uog(e 2my + e~ 2my + e mi? ), 


(4.21) 


D 4 (m) 


: Re^ = 


\og(e my + e- my + + e mR ). 


(4.22) 
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It is interesting to observe that the double wall solution D (m) in Eq.(4.21) in the case of infinite 



mR i n 




8 
6 




4 




(4.26), 




/ (4.24) 



-5 -2.5 2.5 5 7.5 10 



mR 



Figure 3: Relation between R defined in Eq.(4.4) and R. The solid line represents the relation de- 
fined in Eq.(4.26) for the finite coupling solution D A (m) in Eq.(4.22). The dashed line represents 
the relation defined in Eq.(4.24) for the infinite coupling solution D (m) in Eq.(4.21). 



coupling can be rewritten as a superposition of two single walls placed apart by a distance R 

Re^(y) = - log(e mfe -^ + e"" 1 ^! >) + - log(e m ^ + t) + e"" 1 ^)). (4.23) 

This parameter R can be regarded as another choice of a moduli parameter and is related to R 
defined in Eq.(4.4) 



iR i —mR 'mR 



e"" 1 + e~"" 1 = e" 



/e + ve"' m — 4\ 
mi? = log ( j, (4.24) 



which is illustrated in Fig. 3. Unfortunately the new choice of the moduli parameter R becomes 
pure imaginary and loses an intuitive meaning as the distance between the two walls when 
R < log2/m. This situation has some similarity to a moduli parameter for a model of BPS 
double wall in an M = 1 SUSY Wess-Zumino model in four-dimensions [?]. We can also rewrite 
the double wall solution D A {m) in Eq.(4.22) in the case of the finite coupling as a superposition 
of two single walls placed apart by a distance R 

Re^(y) = log(e^-f ) + e -^v-*)) + log(e^+f ) + e -^ {y+ ^). (4.25) 

The new moduli parameter R is related to R in this case as 

e~ +e~— = V6 + e mR mR = 2 log^- ^ J. (4.26) 

In this case of the finite coupling D 4 (m) in Eq.(4.22), the parameter R takes only positive real 
values, whereas R takes positive as well as negative values as illustrated in Fig.3. In both cases of 
infinite and finite coupling, both R and R have an intuitive meaning of relative distance between 
the wall, as long as the distance is large : R — > R, for R — > oo. Therefore R in this case gives an 
intuitively nicer parametrization for the relative distance between the two walls. 
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a) b) 

Figure 4: Comparison of the vector multiplet scalar E as a function of my for exact solutions 
with various gauge couplings, a) Single wall solutions 60 (m) (dotted line), 6 2 (m) (solid line), 
6 3 (ra) (short dashed line), and 64(771) (dashed line), b) Double wall solutions D (m) (dotted 
line), and D^m) (solid line). 

It is interesting to examine how these double wall solutions behave, in the limit of coincident 
walls and in the limit of asymptotically far apart walls. By smoothly changing the moduli 
parameter R, we find the following limiting behaviors which can be expressed symbolically as 

S (m) R ^° D (m) ^S (|)©So(f), (4.27) 
5 4 (m) R ^-°° £4(777) *=±?° 5a (y) © S 2 (y) . (4.28) 

These limiting behaviors apply not only for the ij)(y), but also for all the physical quantities. 
Note that the sum of the numbers k is preserved, corresponding to the fact that the total tension 
is a conserved topological invariant for a given boundary condition. We can find E, H 1A and the 
potential V corresponding to these exact single and double wall solutions (4.17)-(4.22) by using 
(4.1), (4.2) and (3.41). We illustrate and compare various single and double wall solutions for 
scalar E(y) of vector multiplet in Fig.4, for hypermultiplet scalars H 1A (y) for A = 1,2 (single 
wall) and A = 1,2,3 (double wall) in Fig.5 and for potential V(y) in Fig.??. 




-6 -4 -2 2 4 6 -10 -5 5 10 



a) b) 

Figure 5: Comparison of the hypermultiplet scalars H 1A as a function of my for exact solutions 
with various gauge couplings, a) Single wall solutions 60(777) (dotted line), 62(777) (solid line), 
63(777) (short dashed line), and 64(777) (dashed line), b) Double wall solutions D (m) (dotted 
line), and D 4 (m) (solid line). 
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a) b) 
Figure 6: Comparison of potential V = normalized by the asymptotic value (y — > ±oo) as a 
function of my for exact solutions with various gauge couplings, a) Single wall solutions S (m) 
(dotted line), S 2 (m) (solid line), S 3 (m) (short dashed line), and S 4 (m) (dashed line), b) Double 
wall solutions D (m) (dotted line), and _D 4 (m) (solid line). 



4.2 Moduli Dynamics of Two Domain Walls 

The low-energy effective Lagrangian for massless field has been studied systematically, especially 
in the case of infinite gauge coupling [?]. It has been proposed to use an expansion in inverse 
powers of gauge coupling, since the infinite gauge coupling case was the only explicitly known 
solution at that time. We would like to derive the low-energy effective Lagrangian for massless 
field in the case of finite gauge coupling, using our exact solutions. We will compare it to the 
previous result of infinite gauge coupling, and show an instability of an expansion in inverse 
powers of gauge coupling. 

Since the overall position and phase become just free massless fields [?], we shall concentrate on 
the relative distance R a and relative phase 9 a of multiple walls defined in Eq.(4.4). These moduli 
parameters represent positions in y and in internal space of U(l) of multiple walls. We have 
obtained mode functions of massless Nambu-Goldstone fields by differentiating the background 
in terms of these moduli parameters in Sect. 3. 2. Here we are interested in not only the mass 
spectrum but also the entire effective Lagrangian at low energies. To find out the low-energy 
effective Lagrangian, a more systematic method by Manton [?] is useful. In this method, we 
need to promote the parameters of the solution (moduli) to be functions of the world volume 
coordinates x^ 1 , namely four- dimensional fields. Assuming the variation in the world volume 
coordinates x M to be weak, we can obtain all the nonlinear dynamics containing smallest number 
of derivatives (two derivatives). This well-defined procedure allows us to obtain the nonlinear 
kinetic term of the massless fields. Since we are considering the massless Nambu-Goldstone 
fields, there should be no additional potential terms without derivatives. We can keep only 
t = x° dependence of these fields, since the Lorentz invariance in the world volume allows to 
recover the entire dependence on x M [?]. Effective Lagrangian C cS of moduli fields is found to be 
given by integrating over the Lagrangian of the fundamental theory evaluated by the background 
field with the t dependent moduli r a (t) = R a (t) + i9 a (t) [?] 

roc N 

£ cS = / dyY,\H 1A \\^ + m A y -a a A n(rn A yl-a a A fi)- (4.29) 

J-oc A=1 

Let us take two wall case (N = 3) for concreteness. We discard terms for the overall position 
and phase, and omit the subscript a, since there is only one set of relative moduli R, 9 now. We 
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will be interested in low-energy effective Lagrangian for the relative distance R and the relative 
phase 9 of two walls. To obtain the dependence on the relative phase 9, we need to find the 
imaginary part of i[)(y) besides the real part that we have determined. The BPS equation (4.2) 
shows that the imaginary part of ip(y) is given by integrating the extra dimension component 
of the gauge field : d y lmip(y) = W 4 (y). As shown in the case of the single wall in Sect. 3. 2, we 
have, near the walls, a nontrivial field strength F m {W) ^ proportional to the derivative of the 
moduli fields 9. Then the extra dimension component of the vector potential is determined by 
the field equation for the vector potential [?]. The most convenient gauge to make the remaining 
SUSY manifest is W„ = and W 4 ^0. 

We find from the equation of motion that the imaginary part (combined with the real part) 
satisfies 

g-Vc(*-,^)e— . (4.30) 

Because of the boundary condition ijj = at y = ±oo, the solutions of the real and imaginary 
parts are proportional to R and 9 with the same coefficient which is a function of R only 

implying W 4 = d y -J^Reijj(y : R)9. The moduli field r = R + iO is a complex scalar of a chiral 
scalar field, taking values [?] R e R and 9 G [0,47r/m). We see that the M = 1 SUSY in 
four dimensions is now manifest. Moreover, the Kahler metric /C rr * of the low-energy effective 
Lagrangian £ e g is given in terms of F(R) which is a function of R only 

CcK= (i? 2 + ^ 2 )/C rr *(r,r*), (4.32) 
/C rr .(r,r*) = ^-m(F(mR) = ^-m(F ( rn-^— ) . (4.33) 



The function F(R) is most conveniently evaluated by taking 9{t) = and by writing hypermul- 
tiplets H lA in terms of the function ip(y) in Eq.(4.1) and using the choice in Eq.(4.13) 

-m(R 2 F(mR) = / dy ^ \H 1A \ 2 (-Reija A R + a 2 A R 2 ) 

A=l 

1 r\ r°° 

= - A mCR 2 iB / dye~ 2 ^ +mR . (4.34) 



4 ' dR 



oo 



The Kahler metric K. rr * gives equations of motion for the relative distance R and for the 
relative phase 9 

R = _I m (iogF)'(i? 2 - 9 2 ), (4.35) 

9 = - m (\ogF)'R0, (4.36) 

respectively, where ' means d/dz, z = mR. These equations of motion are real and imaginary 
parts of a single equation due to the complex structure of the N = 1 SUSY : 

f = --m(logF) V 2 , r = R + i9 (4.37) 
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We observe that the above equations of motion admit that 9 = is always a consistent solution. 
Therefore we can discuss the relative motion of the double wall with fixed relative phase consis- 
tently. On the other hand, in order to have a fixed R, the above equations of motion is consistent 
only if is also constant. Therefore we shall consider the relative motion assuming a constant 
relative phase 9 in the following. 

Inserting our exact wall solution D^{m) with the finite gauge coupling in Eq.(4.22) into 
Eq.(??), we obtain the corresponding Kahler metric F(z = mR) for the finite gauge coupling 



F(z) 



dt 



e* + e - * + 



(e* + e-* + V6 + e*) 3 
24 



4 + 



log 



^6 + 6^ + ^2 + e 2 



]}■ 



(4.38) 



(2 + e 2 )H v/(2 + e*)(6 + e*) 

For the D (m) solution at the infinite coupling, the Kahler metric F has been obtained as [? 




-10 



10 



Figure 7: Comparison of Kahler metrics F(z) for the infinite gauge coupling (dashed line) and 
F(z) A for the finite gauge coupling (solid line), as functions of z = mR. 



F(z) = - dt 



e* r . e* + e _i 



2 J-oo (e* + e-* + e z ) 2 

log [ 1 ) . (4.39) 



e 22 - 4 I Ve 22 - 4 U 2 + ^ 

The metric is real and positive for the entire values of — oo < z < oo, in spite of the apparent 
singularity at e 2z = 4 [?]. These Kahler metrics are illustrated and compared in Fig.??. We see 
that there is a significant difference between them quantitatively, although general features are 
similar. 

The physical significance of the quantitative difference can perhaps most easily appreciated 
by comparing the strengths of forces acting between the walls in the case of the finite coupling to 
that of the infinite coupling. In Fig.??, we illustrate and compare the coefficients F'/F of force 
in Eqs. (??)-(??) for the infinite and finite coupling cases. We observe that the strength of the 
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force for the finite coupling case is much larger than the infinite coupling case for large relative 
distance R. 



F{z)'/F{z) 




Figure 8: Comparison of the ratio F' Q /F (dashed line) and F'JF± (solid line). 



To understand the quantitative difference of the strength of the forces at large values of R, 
let us evaluate the Kahler metric F at asymptotic values of R. In the region of R — > — oo, both 
target space metrics F(z) are flat and are related each other by a scale transformation. On the 
other hand, in the case of R — > +00, our Lagrangian for the finite coupling becomes 



£4 = -^rnC 



1 - 8e~ mR + (!2mR + 28)e" 2mK {R 2 + 9 2 ). 



2mR 



(4.40) 



Consequently the equation of motion of the relative distance R for D 4(772) solution at the finite 
gauge coupling is found to be 



R 



-e~ mR Am + e- 2mR (12m 2 R - 10m) (R 2 - 9 2 



On the other hand, the asymptotic Lagrangian for infinite coupling at R 

to be m 



1 



1 - A(mR - l)e~ 2mR (R 2 + 9 2 ), 



(4.41) 

00 has been found 
(4.42) 



and the corresponding R — > 00 limit of the equation of motion of the relative distance R for 
D (m) solution at the infinite coupling reads [?] 



R = 



-2mR 



{Am 2 R-Qm){R 2 -9 2 ). 



(4.43) 



We see that for R 2 > 9 2 there is always an attractive force operating between walls which are 
sufficiently far apart, irrespective of the strength of the gauge coupling. However, we find that 
the strength of the force for the finite coupling behaves as e~ mR , which is much larger than that 
for the infinite coupling case of e~ 2mfi . This explains the reason why we obtained much stronger 
force between the walls in the case of the finite coupling. 
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4.3 Approximation in 1/g 2 and Asymptotic Behaviors 



Here we wish to discuss the previously proposed expansion in inverse powers of gauge coupling 
to obtain finite coupling results [?]. We also analyze the asymptotic behavior of the BPS wall 
solutions for generic values of gauge coupling to understand the reason why the power series 
approximation exhibits pathological behavior. 

We expand Re^ in inverse powers of the gauge coupling : a = (2g 2 £) _1 = /c 2 /(4m) 2 as 

oo 

Reifj = a n ^n- (4.44) 

n=0 

By substituting the expansion to the equation (4.5), we obtain an expansion in power of a. The 
leading power 1/a comes only from the right-hand side and gives the result at the infinite coupling 
in Eq.(4.21) : ip = W. Using this result, all the successive powers can be solved iteratively 

00 00 ^2 / oo / n / oo \ 

E«**^E«"^- 2 E^w E«'* + , • (4.45) 

n=0 n=0 y n=l V ; ' \l=0 / 

Several lower order results are given explicitly by 

fa = W, Vi = W {2 \ i> 2 = + {W {2) )\ 

+ \2{W^fW {2) + 4(H/ (2) ) 4 , • • • , (4.46) 

where is defined by = d n W/dy n . 

First let us compare the power series approximation with our exact solution of the single wall 
5-2 (m) for the finite gauge coupling k = 2. We illustrate the result of the power series expansion 
up to the Z-th order S^- ) (y) = d y ip^(y) = X^L=o ^^y^niv) by setting a = (2g 2 Q~ l = (2m) -2 , and 
compare it with the exact solution d y ijj(y) in Fig.??. Near the wall, the approximation oscillates 
wildly and shows no indication of convergence. 



Similarly we can obtain the power series approximations for the Kahler metric F for the 
relative position and phase for two walls. For concreteness, we will take the case of our exact 
solution D±{m) for the finite coupling k — 4, and compare the exact Kahler metric F(z) 4 with 
a profile obtained by the approximation in inverse powers of the gauge coupling. We can obtain 
the approximations up to the Z-th order F(z)f\ by inserting the above expansion (??) to the 
action (??) and by setting a = (^g 2 ^)^ 1 = mr 2 . In Fig.??a), we compare the power series 
approximation up to various orders starting from the zero-th order F (infinite gauge coupling 
case) up to the third order approximations F^ 3 \ We see wild oscillations near the wall. As 
the order of approximation increases, this oscillatory behavior becomes wilder, without any 
indication of dying out. If we abandon to describe the region near the wall, we can examine 
the approximation in the asymptotic regions away from the wall. We also illustrate the Kahler 
metric F(z) for larger values of z — mR — > oo in Fig.??b). Even there, there is no indication 
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Figure 9: Comparison of approximations in inverse powers of gauge coupling to the exact solution 
S ex (|/) for the vector multiplet scalar E(y) = d y Keijj(y) with a finite gauge coupling g 2 ( = 2m 2 . 
Approximation up to the Z-th order is denoted as E^. The zero-th order approximation (infinite 
gauge coupling) E( ) is represented by a dashed line. E^ by a dash-dotted line, E^ 2 ) by a short 
dashed line, and E^ 3 ) by a dotted line. The mass parameter is set to m — 1. 

for the series to converge to our exact result. Actually, we can make the discrepancy more 
quantitative. For instance let us take the power series approximation up to the fourth order 
F(mR)^ for the function F(mR) with an arbitrary coupling a = (2g 2 ()^ 1 = (k/(4m)) 2 . Its 
asymptotic behavior at R — * oo is given by 

F(mR)f ] = 1 + (-4(1 + Am 2 a + 16m V + 64m 6 a 3 + 256m 8 a 4 )mR 

+ 4 + 24m 2 a + 128mV + ^m 6 a 3 + ^^mV + 0(a 5 )) e~ 2mR 

+0(e- 4mR ), (4.47) 

whereas the exact form of the function with the finite coupling at R — > oo is 

F{mR) 4 = 1 - 8e~ rnR + (12mR + 28)e" 2mR + 0{ e - 4mR ). (4.48) 

It is now clear that the power series approximation F(mR)^\l — > oo can never converge to the 
function F(mR)±. Therefore we conclude that the approximation by an expansion in inverse pow- 
ers of the gauge coupling is pathological, and may not be suitable to extract physical quantities 
correctly. 

Finally we analyze the asymptotic behavior of the BPS solution for arbitrary finite gauge 
coupling. This can be worked out even without exact solutions and may provide an alternative 
approximation scheme, rather than the expansion in inverse powers of gauge coupling. Let us 
take single wall case with two hypermultiplets for simplicity. Since BPS solutions are essentially 
real when the moduli parameters are constants (not fields), we can ignore Inr0(y) and obtain the 
BPS equation as 

d 2 ^j 



= g 2 ((l-e- 2 ^{e 2my + e' 2my )) (4.49) 
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F(z) 



0.98 
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b) 



Figure 10: Comparison of approximations in inverse powers of gauge coupling up to the Z-th 
order F±\z) to the exact solution F 4 (z) for the Kahler metric of moduli fields as a function of 
z = mR. a) Near the region of coincident walls, b) Asymptotic region y = mR — > oo. 

Let us parametrize the gauge coupling in terms of k which is no longer an integer here : g 2 ( = 
8m 2 /k 2 . We wish to devise an approximation scheme valid for large values of z = my. Since 
the vector multiplet scalar E = d y ip(y) should approach to a constant asymptotically, the BPS 
equation dictates that ip{y) should approach ±my at y — > ±oo. We expand the remaining 
subleading order terms in a series which will be determined successively as solutions of iterative 
equations 



j oo 

if>(y) = my + - ^2e n {y), 



(4.50) 



n=l 



where we can assume that the remaining terms should vanish asymptotically, since we can absorb 
a possible constant term into a shift of y 



£ n (y)^0 asy^O, £i(y)>£ 2 (y)> 
By substituting the expansion to the BPS equation, we obtain 



d 2 e r 



1 v ^ <• - 

n=l a 



:i + e 



-Amy 



(4.51) 



(4.52) 



We can rewrite this equation with z = my and g 2 ( = 8m 2 /k 2 



E 

n=l 



^(l-e- &£ "(l + e- 4 *)), 

-e- 42 (l- £l + 



16 / 1 2 

£i + ^2 - -^i + 



k 2 ' ~* 2 

The first order approximation is determined by the following equation 

d 2 e x 16 



dz 2 



k 2 



(ei - e 



-42 



). 



(4.53) 



(4.54) 
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which yields a general solution with an arbitrary integration constant C\ 



Az 1 _A-y 

£l = 2 Cl e * + Y^e~ Az , 



(4.55) 



where we have used the boundary condition e± — > at y — > oo to eliminate another possible 
integration constant. Therefore the asymptotic behavior of the first order approximation Si{y) 
depends on the strength of the gauge coupling constant g 2 = j^- as 



The wall solution requires that the vacuum value m is reached from below : = d y ifj(y) < m 
as y — > oo, This implies £i(y) > 0. Therefore we obtain up to the first order approximation 



This result can be regarded as the origin of a different asymptotic behavior for the finite coupling 
exact solution (k — 2 case) compared to the infinite coupling solution (k — case). It is 
interesting to observe that the critical case of k — 1 has an exceptional asymptotic behavior 
corresponding to degenerate exponents. This complication may be related to the fact that we 
have not yet succeeded to obtain an exact solution in that case. 

It is now straightforward to extend our approximation scheme to more general multi-wall 
cases. Since our approximation scheme from asymptotic region is applicable to any values of 
gauge coupling constant, this should be useful as an alternative approximation method instead 
of the expansion in inverse powers of gauge coupling. Our iterative approximation is at least 
guaranteed to be valid for asymptotic region y — > oo, although it may wildly oscillate and possibly 
diverge near the wall. On the other hand, we have seen that the expansion in power of inverse 
gauge coupling does not converge near the wall nor asymptotically. To obtain a full solution, one 
needs to determine the integration constant such as c\ which can be done by smoothly connecting 
to solutions from the region near the wall. Only when this integration constant is chosen to be 
a particular value, the solution should smoothly connect to the solution near the wall. 
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